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ABSTRACT 

Since  a  nonlinear  feedback  control  system  may  possess  more  than  one 
type  of  forced  oscillations,  it  is  highly  desirable  to  investigate  the 
type  of  forced  oscillations  which  can  occur  when  the  nonlinear  restoring 
force  function  is  of  a  specific  form. 

This  paper  is  to  present  a  harmonic  linearization  method  for  finding 
the  existence  of  forced  oscillations  and  response  curve  characteristics 
of  a  nonlinear  feedback  control  system  by  means  of  finding  the  restoring 
force  function  of  the  nonlinearity  and  using  the  harmonic  balance  and  an 
iteration  method  for  investigating  the  conditions  for  one  type  of  forced 
oscillations  exhibited. 

The  existence  conditions  for  fundamental  frequency,  sub-harmonic  of 
2nd  order  and  3rd  order  forced  oscillations  of  a  second  order  feedback 
control  system  are  investigated;  also  the  fundamental  frequency  forced 
oscillation  for  a  higher  order  system  and  the  jump  resonance  frequencies 
of  response  curve  are  investigated,  and  a  general  expression  for  the  equi- 
valent gain  of  the  nonlinearity  has  been  developed. 

The  author  wishes  to  express  his  appreciation  for  the  assistance  and 
encouragement  given  by  Dr.  George  Julius  Thaler  of  the  U.  S.  Naval  Post- 
graduate School  in  this  investigation. 
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CHAPTER  I 
INTRODUCTION 
One  of  the  major  branches  of  servo  mechanism  analysis  for  nonlinear 
control  systems  is  the  study  and  methods  of  prevention  of  continuous  oscil- 
lation.  For  a  linear  control  system,  the  output  may  oscillate  with  con- 
tinuous or  increasing  amplitude  with  no  input  signal  applied,  in  which  case, 
the  system  is  said  to  be  unstable,  or  they  may  have  oscillations  which  die 
away.   In  both  cases  this  transient  response  is  the  characteristic  of  the 
system  itself  and  the  conditions  for  divergent,  continuous  or  decaying 
oscillation  do  not  depend  on  the  form  or  magnitude  of  the  input  signal. 

When  a  nonlinear  element  is  present,  this  independence  of  the  input 
signal  r.o  longer  holds.   The  characteristic  performance  of  the  nonlinear 
control  system  will  depend  on  both,  i.e.,  forcing  function  and  the 
characteristic  of  the  nonlinear  element. 

A  nonlinear  control  system,  under  suitable  conditions  may  exhibit 
steady  oscillations  in  which  the  main  component  has  a  frequency  which  is 
dependent  on  the  frequency  of  the  forcing  function.   This  type  of  oscilla- 
tion is  called  "Forced  Oscillation".   A  forced  oscillation  for  which  the 
frequency  is  a  fraction  of  the  forcing  function  frequency  is  called  "Sub- 
harmonic  Forced  Oscillation".   A  forced  oscillation  for  which  the  freqtiency 
is  multiple  of  the  forcing  frequency  is  called  "Super  Harmonic  Forced  Oscil- 
lation". 

Since  a  nonlinear  control  system  may  possess  more  than  one  type  of 
forced  oscillations,  it  is  highly  desirable  to  investigate  the  type  of 
forced  oscillations  which  can  occur  when  the  nonlinear  restoring  force 
function  is  of  a  specific  form. 

Some  methods  of  investigation  for  nonlinear  forced  oscillation  have 


been  discussed  by  J.  C.  West   tn  his  book  "Analytical  Techniques  for 
Nonlinear  Control  System"  (1960).  A  graphic  method  for  investigating 

forced  oscillation  of  "ON  OFF"  control  system  has  been  developed  by  Hamel 

2 
(French)  and  Tsypkin  (USSR)  .   Some  work  for  subharmonic  oscillations  of 

a  typical  restoring  force  function  has  been  done  by  C.  A.  Ludeke  and 

William  Pong  in  their  paper  1959,  and  also  by  Ogata  in  his  PhD  thesis  , 

June  1956. 

There  are  some  other  authors  for  investigation  of  the  vibration  of  a 

5  6  7 

mechanical  system,  such  as  N.  Minorsky  ,  J.  J.  Stoker  ,  and  Y.  H.  Ku  in 

their  books  of  "Introduction  to  Nonlinear  Mechanics" ;"Nonlinear  Vibrations" 

and  "Analysis  and  Control  of  Nonlinear  System"  respectively  => 

The  purpose  of  this  paper  is  to  present  a  harmonic  linearization 

Q 

method  for  finding  the  existence  of  forced  oscillations  and  response 

curve  characteristics  of  a  nonlinear  feedback  control  system,  by  means  of 

9  10 
finding  the  restoring  force  function  of  the  nonlinearity  9   ,  and  using  an 

iteration  method  for  investigating  the  condition  for  one  type  of  forced 

oscillations  exhibited.  A  response  curve  and  "Jump  Resonance"  can  be 

also  investigated. 

The  emphasis  of  this  paper  is  placed  on  showing  the  general  approach 
of  the  method,  therefore,  second  order  nonlinear  feedback  control  systems 
are  used  as  examples.  However,  a  higher  order  system  will  be  discussed  in 
this  paper  also. 

By  the  harmonic  linearization  and  iteration  method,  a  general  equa- 
tion for  the  equivalent  gain  of  an  odd  function  nonlinear  element  will  be 
developed  .   It  is  in  terms  of  the  amplitude  and  the  frequency  of  input. 

There  are  six  chapters  in  the  main  body  of  this  paper.  The  first 
chapter  and  the  last  chapter  are  a  general  introduction  and  conclusion 


respectively.   Chapter  two  is  a  general  description  of  oscillation  of 
feedback  control  systems  and  an  investigation  of  the  restoring  force  func- 
tion of  the  nonlinear  element.   Chapter  three  investigates  the  existence  of 
fundamental  forced  oscillations  and  the  jump  phenomena  in  the  nonlinear 
feedback  control  system.   Chapter  four  is  a  sub  harmonic  oscillation  in- 
vestigation.  Chapter  five  investigates  the  existence  of  forced  oscillation 
of  a  higher  order  feedback  control  system. 


CHAPTER  II 

GENERAL  DISCUSSION  OF  OSCILLATION  IN  FEEDBACK  CONTROL  SYSTEMS 
2-1  General  Description: 

For  either  a  linear  or  nonlinear  feedback  control  system,,  most  re- 
quirements for  control  system  design  specify  a  minimum  time  response,  high 
accuracy,,  and  high  stability.   It  is  clearly  in  order  to  have  a  minimum 
time  response  and  high  accuracy  for  a  feedback  control  system,  but  first 
of  all,  the  system  should  be  stable. 

In  general,  two  types  of  oscillations  may  be  exhibited  in  a  feedback 
control  system,  one  is  free  oscillation,  i.e.,  an  autonomous  system.   This 
is  a  system  which  is  left  by  itself  with  its  outside  source  of  energy  sudden- 
ly cut  off  or  abruptly  changed  in  amplitude.   In  such  a  system,  free  oscilla- 
tions occur  as  the  system  tends  to  equilibrium  after  reaching  a  state  of  un- 
equilibrium.   The  oscillation  frequency  of  this  type  depends  only  on  the 
characteristic  of  system  itself.   Ordinarily,  a  free  oscillation  is  a  damp- 
ed oscillation;  as  the  time  increases  the  amplitude  decreases  and  damps  out 
in  steady  state.   For  an  ideal  case,  if  it  is  a  system  without  damping,  i.e., 
oC   s>  0,  the  oscillation  will  be  continuous. 

The  other  type  of  oscillation  is  a  forced  oscillation,  i.e.  a  nonauto- 
nomous  system.  This  is  a  system  that  is  acted  upon  by  an  outside  source  of 
energy,  a  force  function.  The  force  function  may  be  a  constant;  a  periodic 
function  or  any  other  function  of  time. 

The  analysis  and  methods  for  preventing  oscillations  of  a  feedback 

control  system  are  concerned  directly  or  indirectly  with  the  characteristic 

of  system,  i.e.,  the  differential  equations  of  the  system.   If  the  system 

is  linear  a  frequency  response  of  the  system  can  be  directly  solved  from 

0 
the  differential  equation,  furthermore,  the  principle  of  superposition  can 

be  applied  too.   On  the  other  hand,  if  the  system  is  nonlinear8  either  with 


a  nonlinear  damping  or  nonlinear  restoring  force  function,,  then  the  basic 
tools  for  the  analysis  of  a  linear  system  are  no  longer  valid. 
2-2  Oscillation  in  Linear  Systems: 

There  are  two  types  of  methods  for  analysis  and  design  of  linear  feed- 
back control  systems;  one  is  graphic  method,  the  other  is  an  analytical 
method.   The  best  known  of  graphic  methods  consist  of  Bode  diagram;  Polar 
plot;  Nichols  diagram;  and  Root  locus.   If  we  know  the  transfer  function  of 
the  system,  any  of  the  above  methods  can  be  used  for  solving  the  problem  of 
stability.   The  analytical  method  is  a  mathematical  analysis  for  solving  the 
differential  equation,  and  plot  the  curves  of  the  solved  equation. 

First  consider  a  simple  linear  feedback  control  system^  as  shown  in 
Fig.  2-1.   The  forward  element  has  a  transfer  function  G(s)s  and  a  feedback 
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Fig.  2-1  Block  Diagram  of  a  Linear  System 
element  transfer  function  H(s).   The  feedback  signal  is  H(s)9c.  (s)  and  it 
is  assumed  that  this  is  subtracted  from  the  input  signal  Qt    (s)  to  form 
the  control  signal  e(s);  thus: 

e(s)  -   9,  (s)  -  H(s)  9c(s)  (2-1) 

The  relation  between  the  control  signal  and  reference  signal  in  the 
Laplace  transformation  form  is: 


a  (s) 


-M*l 


1  +  H(s)G(s) 


t,(s) 


(2-2) 


The  solution  of  equation  (2-2)  for  a  particular  input  can  be  obtained 
by  Laplace  and  Heaviside  methods.  This  involves  the  determination  of  the 
roots  of  equation: 

1  +  H(s)G(s)  =  0  (2-3) 

In  order  to  split  the  transfer  function  into  the  sum  of  partial  fractions 
of  standard  form  for  vhich  time  solutions  are  known.   If  equation  (2-3)  is 
an  n  th  order  polynominal  in  s,  there  will  be  n  roots  say  S  to  S  ,  such 
that  (s-  S  )  is  a  factor  of  the  polynominal.   Special  cases  arise  in  which 
some  of  roots  are  identical,  and  if  m  roots  of  S  occurs,  (s  -  S  )   is  a 
factor  of  the  polynominal. 

In  general,  the  roots  S  will  be  a  complex  number^  as: 

S  ■  ^       +   1  CO  (2-4) 

r      r    J  r 

in  which  ^      and  u1     are  both  real  numbers,  and  since  the  coefficients  of 
r       r 

the  polynominal  of  equation  (2-3)  are  all  necessarily  real,  then  the  com- 
plex roots  should  be  in  conjugate  pairs,  Thus: 

(s  -^r  +  j^r)   (s-«r  -  j  A)  r) 
will  be  a  factor  of  the  polynominal. 

In  determining  the  transient  response  to  a  forcing  function  input, 
these  various  factors  become  the  denominator  of  the  partial  fraction  re- 
presentation of  the  transfer  function  equation  (2-2).  The  term  of  K  /(s  -S  ) 
gives  a  time  solution  containing  expoential  term: 

(^  r  +  j  Co    )t 
ev      J  x' 

The  complete  solution  is  the  sum  of  all  such  terms. 

The  complex  conjugate  pairs  give  rise  to  an  oscillatory  term: 

e  ^Crt       Sin(<iy  t  +  f    ) 
If  the  system  has  convergent  response,  all  of  these  terms  in  the  time  solu- 
tion must  converge  and  hence  all  roots  S.  to  S  must  have  negative  real 
parts.   If  any  of  these  roots  has  a  zero  real  part,  then  a  continuous 


oscillation  is  produced.   If  any  of  the  complex  roots  has  a  positive  real 
part,  then  divergent  oscillation  will  be  the  result. 

To  determine  the  stability  of  a  system  by  using  the  graphic  method,  it 
is  not  necessary  to  find  the  roots  of  the  characteristic  equation.   If  the 
equation  is  known  algebrically,  the  Routh  method  can  be  applied  and  also 
by  the  Nyquist  and  Bode  diagram  method  a  graphical  solution  can  be  obtained. 

Nyquist  criteron  graphical  analysis  is  the  most  useful  method  to  solve 
the  stability  problem  of  control  systems.   Rewriting  equation  (2-3)  in  the 
form  of; 

G(s)H(s)  — 1  (2-5) 

The  roots  of  equation  of  (2-3)  are  particular  values  of  s  which  satisfy 
equation  (2-5).   Thus  if  the  equation  of  G(s)H(s)  can  be  mapped  for  continu- 
ously variable  values  of  s,  then  those  values  for  which  G(s)H(s)  ■  =1  can 
be  determined. 

Assume  a  particular  roots  of  S  ■  U       +  j    can  be  represented  by  a 
particular  value  of  G(s)H(s),  i.e. 

GU  x   +  Kr)H(.ir  +  j   r)  -  R(   ,  :..)   +  jl(   ,    )     (2-6) 

This  equation  represents  a  vector  in  I  vs  R  plane,  the  magnitude  of  vector; 

I  2    2 
M  =  J  R  +  I  (2-7) 


and  the  direction  by  an  angle 


Tan_1l/R  (2-8) 


For  a  constants  ,  the  value  of  M  will  be  a  locus  in  I-R  plane  as 
r  v 

to   varies.   This  process  can  be  repeated  for  different  values  of  '  r  to 
obtain  a  family  of  curves  as  shown  in  Fig.  2-2. 

This  steady  state  frequency  response  locus  of  the  open  loop  system 
can  be  obtained  experimentally  from  the  system  by  measuring  gain  and  phase 
for  a  sinusodial  input  of  varying  frequency.   It  can  be  shown  that  is  the 
steady  state  frequency  response  locus  passes  to  the  right  side  of  (-1,  0) 
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Fig.  2-2  Locus  of  an  Open  Loop  Transfer 
Function  G(s)H(s) 

point,  this  system  will  be  stable,  on  the  other  hand8  if  it  passes 
through  or  to  the  left  side  of  (-1,  0)  point,  the  system  will  be  unstable. 
This  criterion  is  due  to  Nyquist.   A  typical  curve  is  shown  in  Fig.  2-3. 

In  the  extension  to  nonlinear  control  systems,  or  more  complex  multi- 
loop systems,  the  single  loop  criterion  is  sufficient  for  the  majority  of 
needs  where  G(s)  can  represent  the  over  all  transfer  function  of  several 
loops  provided  that  G(s)  is  itself  stable. 


-£2± 


H,e) 


Fig.  2-3  Nyquist  Criterion  for  Three 
Characteristic  States 


2-3  Forced  Oscillations  in  a  Linear  System: 

Assume  a  second  order  linear  system  block  diagram  is  shown  in  Fig. 
2-1;  in  which: 

and: 

H(s)  m   1  (2-10) 

From  equation  (2-2),  the  differential  equation  of  system  becomes: 

Q     +  <S  9       +  KG  «  KG  (2-11) 

c      c       c     r 


Let  the  forcing  function: 


Equation  (2-11)  becomes 


0  (t)  »  F  Coswt  (2-12) 


e   +  /   6   +  KG  *  C  Coswt  (2-13) 

c        c       c 

in  which  C  is  a  constant  of  value  of  FK. 

The  solution  of  equation  (2-13)  consists  of  the  sum  of  the  solution 
of  the  homogeneous  equation^  (i„e.  the  free  oscillation  of  the  system)  and 
the  solution  of  the  non-homogeneous  equation.   Let  the  solution  of  homo- 
geneous equation: 

9c(t)  «  e"*  fc(C,  CoscOn  J   l-/t  +  C2  Sin  WnJ    l-f  t)   (2-14) 
in  which: 

60  n  -  JY 
J.     «  2  U'.t. 
f     *  Cos  I' 
The  complete  solution  of  equation  (2-13)  will  be: 

0  (t)  *  e  ~"> fc(C.  Cos    Jl^r   t  +  C,  Sin     1  ;'t)  *   JL^^ t _+_ Q)  (2-15) 

C  in  l  nv        ^  (rf-ui'Ti- +/%*&■ 

Equation  (2=15)  is  obtained  by  a  superposition  of  the  free  oscillation 
and  the  forced  oscillation  which  varies  from  the  action  of  external  force: 

The  frequency  of  forced  oscillation  is  the  same  as  that  of  the  external 


force,  the  magnitude  of  the  forced  oscillation  is  given  by: 

H  -  — (2-16) 

The  value  of  phase  shift  0  related  to  the  external  force  is  given 
by: 

CosG  -        tV*  '  ^  (2-17) 

■  ■ 

SinG  ^  (2-18) 

In  the  case  of  positive  damping,  i.e.  U  y   0,  it  is  clear  from  equation 
(2-15)  that  after  a  sufficiently  long  period  of  time,  the  free  oscillation 
is  damped  out  and  only  the  forced  oscillation  would  be  observed. 

In  the  case  of  no  damping,  i.e.  a    ■  0,  the  phase  shift  0   is  seen 
from  equation  (2-17)  and  (2-18)  to  be  zero  for  w<a)n   and  tl   for  cU^uJrf      In 
other  words,  the  forced  oscillation  is  in  phase  with  the  external  force  if 
the  forced  oscillation  frequency  is  less  than  the  free  oscillation  fre- 
quency, and  is  180  degrees  out  of  phase  with  the  external  force  when  u)   is 
greater  than  tCn  . 

In  the  case  of  X   »  0,  equation  (2-15)  becomes: 

6  (t)  »  (C,  Cos  ^  t  +  C0  Sint/;  t) 
c        1      n     2      n 

+   C  Cos(  -  t  +  9)  (2-19) 

I '      !  I 

There  will  be  a  superposition  oscillation  of  two  frequencies,  if  the  value 

of  co  %-  o>*j  one  of  which  is  the  natural  frequency,  and  the  other  is  the 

frequency  of  external  force.   In  the  case  of  a)     *„ ,  ,  the  free  and  the  forced 

oscillations  have  the  same  frequency,  and  the  solution  of  equation  (2-13) 

will  be  found: 

C 
6c(t)  »  C  Cosh  t  +  C2  Sin.  t  -  -^-t  Sintvt  (2-20) 

In  this  case,  the  component  due  to  the  external  force  is  no  longer  periodic, 
it  is  oscillatory  with  an  amplitude  that  increases  linearly  with  time.   At 
this  condition,  it  is  a  resonance  phenomenon.   It  is  very  important  in  de- 
sign to  avoid  this  phenomenon  in  a  system  with  a  periodic  force  function. 
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In  the  case  of  a  damped  system,  i.e.,,     0,  from  equation  (2-16)  for 

the  steady  state,  the  amplitude  of  forced  oscillation  is  always  finite.   It 

is  possible  to  give9  by  the  use  of  dimensionless  variable,  a  more  general 

significance  to  the  expression  for  the  amplitude  of  forced  oscillation,  let 

the  amplitude: 

C 
|H|  -  M(   R  )  (2-21) 

Therefore  equation  (2-16)  becomes: 

M-  f= = (2-22) 

Where  M  is  defined  to  be  the  magnification  factor.   The  extreme  values  for 

2  2  2 

M  are  attained  for  vO  -  0  and  (   /    )   «  1  -  2  f    .   If  1  -  2  f  <  0, 

2 
there  is  a  maximum  for  a   •  0;  if   1  -  2  J     >   0,  and  tc>0,  there  is  a  maxi- 
mum for  \A)/'a)n  —  aH  ~*fx  and  a  minimum  for  co  m   0»   For  small  values  of  damp- 
ing coefficient,  the  frequency  which  produced  the  maximum  amplitude  is  very 
nearly  the  natural  frequency  of  the  system.   Fig.  2-4  shows  the  forced  oscil- 
lation response  curves  of  a  linear  system  as  a  function  of  t-o/  u)^    with 
various  values  of  f     . 
2-4  Oscillations  in  Nonlinear  Feedback  Control  System: 

As  mentioned  in  section  2-1  of  this  chapter,  when  a  nonlinearity  is 
present  in  a  feedback  control  system,  the  properties  of  proportionality 
and  superposition  are  not  valid.   In  the  linear  system,  these  properties 
involve  the  existence  of  a  transfer  function  and  of  characteristic  fre- 
quencies proper  to  the  system  (frequencies  at  which  the  system  tends  to 
oscillate  with  an  amplitude  depending  on  initial  conditions);  in  the  case 
of  nonlinear  systems,  the  amplitude,  like  frequency 8  can  depend  at  one  and 
the  same  time  on  both  the  initial  conditions  and  the  system  itself.   Some- 
times, as  in  the  case  of  limit  cycle  both  frequency  and  amplitude  of  input 
are  characteristic  of  the  system  and  independent  of  of  initial  conditions. 
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Hg.2-4     Response  Curve  Tor  a  Linear 


In  many  cases  the  analysis  and  design  of  such  systems  using  linear 
theory  may  produce  an  excellent  result.   This  general  procedure  is  a 
type  of  linearization  procedure,  i.e.,  the  actual  nonlinear  system  is  re- 
placed by  a  linear  system,  which  approximates  it,  and  the  analysis  and  de- 
sign techniques  are  applied  to  the  linear  equivalent.   The  most  useful  method 
for  nonlinear  linearization  is  that  of  the  "Describing  Function"  method.   It 
is  actually  a  first  harmonic  approximation  method. 

The  method  for  linearization  used  by  this  paper  is  a  multi-harmonic 
linearization  and  iteration  method;  first  assume  a  harmonic  solution  for 
the  system  differential  equation  and  insert  it  in  the  differential  equation, 
and  compare  the  coefficients  of  the  same  order  harmonic  terms „   A  detailed 
procedure  will  be  discussed  in  the  next  chapter. 

As  previously  explained,  a  nonlinear  feedback  control  system  may  posses 
some  type  of  oscillation.   Which  type  will  exist  depends  on  the  restoring 
force  function  of  nonlinear  element  and  the  amplitude  of  forcing  function. 
Hence,  before  investigating  the  existence  of  forced  oscillations,  the  re- 
storing forcing  function  of  the  nonlinearity  should  be  investigated. 
2-5  Characteristic  of  Nonlinear  Elements: 

Nonlinear  elements  which  may  be  seen  in  feedback  control  systems  can  be 
represented  by  a  combination  of  four  fundamental  concepts  of  nonlinearity, 
which  are:  (1)  Saturation;  (2)  Variable  gain;  (3)  Dead  zone;  (4)  Hysteresis. 
These  basic  notions  of  nonlinearity  can,  in  combination  with  one  another, 
result  in  almost  all  types  of  nonlinearity  in  practice,  as  shown  in  Fig.  2-5. 
For  example  an  ideal  relay  nonlinearity  can  be  seen  as  a  saturation  nonline- 
arity which  linear  part  with  a  value  of  infinite  slope. 

The  operational  characteristic  of  a  nonlinear  element  which  operates  in 
a  control  system  depends  on  the  operating  condition9  just  as  a  vacuum  tube 
for  which  the  operating  characteristic  depends  on  the  operating  point  of  its 
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characteristic  curves.   As  a  saturation  nonlinearity  element  with  a  small 
forcing  function  input,  it  operates  linearly.   On  the  other  hand,  if  the 
forcing  function  is  very  large  and  with  a  high  frequency,  the  actual  opera- 
ting characteristic  is  like  an  ideal  relay.   Hence,  the  actual  operating 
characteristic  of  a  nonlinear  element  may  be  said  to  be  a  function  of  the 
amplitude  and  frequency  of  the  forcing  function. 
2-6   Restoring  Force  Function  Investigation  of  a  Nonlinear  Element: 

Consider  the  block  diagram  of  Fig.  2-6,  in  which  the  output  is  a  func- 
tion of  input,  usually  this  function  in  a  control  system  is  called  the 


Nonlinear 


Element 


i 


Fig.  2-6  Basic  Relation  Between  Input 
and  Output  of  a  Nonlinearty 


"Restoring  Force  Function". 
Assume  the  input  is: 


e.  -  E(t) 


The  output  will  be  a  function  of  input,  as 


e     «  f(E) 
o 


(2-23) 


(2-24) 


From  the  numerical  analysis  and  curve  fitting  process,  the  function 
of  f(E)  can  be  expressed: 


f(E)  «  a  +  a.E  +  a0E  +  a,E3+ 
o     1      d.  3 


.+  a 


n 


+  bjE  +  b2E1/2  +  b  E1/2  +  ••••••••••••••••+  bnE 
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1/n 


(2=25) 


Equation  (2-25)  can  be  in  the  general  form: 


f(E)  -  a  + 
7    o 


I  a.El  +     °    b.E1/l 

i  »  1    1        1-1 


(2-25a) 


for  this  paper,  the  nonlinear  characteristic  is  considered  symmetrical  to 
the  original  point,  that  means: 

f(-E)  -  -f(E)  (2-26) 

It  is  an  odd  function,  the  coefficients  of  even  order  of  equation  (2-25)  are 
zero,  it  becomes: 


3      5 
f(E)  -  a  +  a,E  +  a.E  +  a_E  + 
o    1     3      5 


+  b  E  +  b3E1/3  +  b  E1/5  + 


•  o  «  e  e  •  • 


(2-27) 


For  the  saturation  and  variable  gain  nonlinearity,  J.  C.  West  uses 
the  restoring  force  function  as  a  form: 

f(E)  =  a  E  +  a3E3  (2-28) 

in  which  a.  is  greater  or  equal  to  zero,  and  a   is  either  greater^  or  less, 
or  equal  to  zero  corresponding  to  the  case  of  a  "Hard  Spring",  "Soft  Spring", 
and  "Linear  Spring"  saturation  nonlinearity  characteristic  respectively. 
A  set  curves  for  all  values  of  a   is  shown  in  Fig.  2-7. 


0 


Fig.  2-7  Characteristic  curves  of  Eq.  (2-28) 
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In  order  to  generalize  for  the  condition  of  saturation  nonlinearity 
with  dead  zone,  equation  (2-28)  can  be  written  with  a  constant  as: 

.3 


f(E)  «  -a  Sign  E  +  a.E  +  a.E' 
o         l     J 


(2-29) 

if  the  dead  zone  is  less  than  one,  the  characteristic  curve  can  be  shown 
in  Fig.  2-8. 


/ 
/ 
v 

-&o 


Fig.  2-8  Characteristic  Curves  of  Eq.  (2-29). 
As  discussed  previously,  when  a  saturation  nonlinearity  is  operating 
with  a  larger  amplitude  than  the  saturation  voltage  and  a  high  frequency, 
the  operating  characteristic  is  like  an  ideal  relay.   On  the  other  hand, 
under  some  suitable  conditions  a  relay  nonlinearity  control  system  may  ex- 
hibit some  sub-harmonic  forced  oscillations,  as  proved  by  the  paper  of  A.  M. 

12 
Hopkin  and  K.  Ogata   .   This  result  is  the  same  as  a  saturation  nonlinearity, 

Actually,  the  operation  of  a  relay  is  not  a  perfect  discontinuous 
characteristic  element,  when  the  relay  is  operating  from  static  to  closed, 
there  will  be  a  little  amount  of  time  delay  from  starting  to  close  con- 
tacts.  Therefore,  the  operating  characteristic  of  an  ideal  relay  may  be 
considered  a  saturation  nonlinearity  with  a  very  large  slope  of  the  linear 
part,  a  sketch  of  it  is  shown  in  Fig.  2-9. 
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Fig.  2-9  Operating  Characteristic  of 
an  Ideal  Relay 

From  the  above  discussion,  and  for  a  proper  value  of  constants,  the 

restoring  forcing  function  can  be  expressed  as: 

(2-30) 

Similarly,   if  we  consider  the  dead  zone  the  equation  of  (2-30)  can 

be  expressed: 

(2-31) 


f(E)  =  axE  +  a3E3  +  a5E5 


3      5 
f(E)  ■  -a   SignE  +  a  E  +  a  E  +  a  E 


13 


M.  J.  Abzug  *  investigates  the  restoring  forcing  function  of  an  ideal 


relay  nonlinearity  by  a  fifth  root  or  cube  root  of  the  input  as; 

.1/5 


f(E)  -  b5E 


(2-32) 


or; 


f(E)  -  b3E 


1/3 


(2-32) 

A  more  detailed  discussion  about  frequency  response  by  using  the  re- 
storing forcing  function  will  be  covered  in  the  next  two  chapters. 
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CHAPTER  III 
FUNDAMENTAL  FREQUENCY  FORCED  OSCILLATION  IN  NONLINEAR  CONTROL  SYSTEMS 
3-1  General  Description: 

Consider  a  system  shown  in  Fig.  3-1,  in  which  G.(s)  is  the  transfer 
function  of  a  controller  element.   G~(s)  is  the  transfer  function  of  a 
motor  and  gear  element,  G  is  a  nonlinear  element.   H(s)  denotes  a  transfer 
function  of  a  feedback  network. 

When  a  periodic  forcing  function  is  applied  to  the  input,  the  output 
of  the  system  may  be  or  may  not  be  a  periodic  function.   If  it  is  a  peri- 
odic output  and  the  frequency  is  the  same  as  the  frequency  of  the  input 
forcing  function,  it  is  said  that  the  system  is  operating  in  forced  oscil- 
lation with  the  fundamental  frequency. 


/Ov  B>) 


T 


$,(s) 

1 1 

H<*) 


Fig.  3-1  Block  Diagram  of  a  Feedback 
Control  System 

There  are  some  methods  to  represent  the  forced  oscillation  in  a  non- 
linear system.   If  we  keep  an  input  forcing  function  with  constant  frequency 
and  vary  the  amplitude  of  input,  a  response  curve  can  be  shown  in  the  "In- 
put vs  Error"  plane,  Fig.  3-2,  in  which  there  is  a  cut  off  amplitude;  that 

means  when  the  input  is  smaller  than  0   ,  there  is  no  forced  oscillation. 

r-m 

The  other  representation  is  with  the  input  amplitude  constant  vary 
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Fig.  3-2   Forced  Oscillation  Represented 

in  the  E  vs  0  Plane 
r 

the  frequency  of  input.   It  is  shown  in  "Error  vs  Frequency"  plane.   Fig. 
3-2a  is  an  example,  it  is  actually  a  closed  loop  frequency  response  of  the 
system. 


CO 


Fig.  3-2a  Forced  Oscillation  Represented 
in  the  E  vs  iO      Plane 

Either  represented  in  E  vs  ©  or  E  vs  CO  plane,  both  of  them  can  be 

represented  in  a  phase  plane  (E  vs  E) .   A  typical  phase  plane  of  forced 

oscillation  from  the  computer  is  shown  in  Fig.  3-3,  in  which  it  is  represent' 

ed  by  an  ellipse.   The  size  of  it  is  changed  with  the  frequency. 

3-2  Basic  Equation  for  Forced  Oscillation  of  2nd  Order  Nonlinear  Feedback 

Control  System: 

Consider  a  2nd  order  feedback  control  system  shown  in  Fig.  3-4: 
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Fig.   3-3  Forced  Oscil- 
lation represented 
m  Phase  Plane. 
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Fig.  3-4  Block  Diagram  of  a  2nd  Order 
Feedback  Control  System 


The  output  equation  in  Laplace  form: 


Vs) 


Kf(E) 


s(S+^  ) 


in  the  differential  form: 


0  +  .  0   -  Kf(E) 
c      c 


(3-1) 


(3-2) 


Where: 


Equation  (3-2)  becomes: 


ec(s)  -  er(s)  -  e(s) 


e  +  i  e  +  Kf(E)  ■  e  +  <  e 

/    r       r 


Assume  the  input  forcing  function: 


0  (t)  «  F  Cos(u.  t  +  9) 


(3-3) 


(3-4) 


(3-5) 


in  which  0  is  the  phase  difference  between  the  error  signal  and  forcing 
function  of  input.   Inserting  equation  (3-5)  in  to  equation  (3-4): 


or: 


Where: 


E  +^E  +  Kf(E)   -  -A  Cos(u,  t  +  6)  — B  Sin(.  t  +  0)  (3-6) 
E  +  act     +  Kf(E)   -  C  Cos(^t  +  *  +  0  +  0)     (3-6a) 


A  *  CO  F 
B  *  aL  a)   F 


/a2  ; 


FriJoL2--*  «J- 


0  ■  Tan"  B/A  =>  Tan  Ufa. 


(3-7) 
(3-8) 
(3-9) 

(3-10) 
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Equation  (3-6)  or  (3-6a)  is  a  basic  differential  equation  of  a  2nd  order 
nonlinear  feedback  control  system,  in  which  f(E)  is  the  nonlinear  restor- 
ing force  function,  and  K  is  a  constant  of  system. 
3-3  Forced  Oscillation;  Saturation  Nonlinearity  Without  Damping: 

Consider  a  system  without  damping,  i.e.  i  ■  0,  equation  (3-6)  becomes: 
E  +  Kf(E)  -  -A  Cos(  <  t  +  0)  (3-11) 

Let  the  restoring  force  function: 


f(E)  =  axE  +  a3E3  (2-28) 


Equation  <3-ll) : 


E  +  K(axE  +  a  E3)  -  -A  Cos(a  t  +  e)  (3-12) 

7  8 
By  the  harmonic  iteration  method   '  ,  the  solution  of  equation  (3-12)  can 

be  expressed  as  a  series  of  odd  order  periodic  harmonics: 

e(t)  -  E.  CosoJt  +  E  Cos3  t  +  E  Costa  t (3-13) 

in  which  E..,  E  ,  E_  are  the  amplitude  of  harmonics  of  e(t). 

If  we  consider  only  the  fundamental  frequency  forced  oscillation,  we 

can  assume  one  solution  is: 

e(t)  »  E-  Cosat  (3-14) 

Where  E.  is  the  fundamental  frequency  amplitude  of  error  signal,  CO   is  the 

frequency  of  forcing  function  and  forced  oscillation;  is  to  be  determined,, 

Since: 

E     »  -   uTE     Cos  a:  t  (3-15) 

E     «  -    bl    E     Sin.Ot  (3-16) 

E3  «  £  E^O  Cos^t  +  Cos3 fl)t)  (3-17) 

Inserting  equation  (3-14)  into  (3-17)  in  equation  (3-12): 

2  3      3  13 

(- a,  Ex  +  Ka1EJ  +  £  Ka3Ej  )  Cos  U)t   +  *  a  E  Cos3^t  (3-18) 

■  -A  (CosiUt  CosG  +  Sin^'t   Sin9) 
neglecting  the  higher  order  harmonic  terms  and  equating  the  coefficients 
of  Cos^t  and  SinU/t,  thus: 
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(-tX>  E,    +  Ka,E,   +7     Ka^El3)    -   -  ACos©  (3-19) 

1114  3 

ASln  n>t      Sine     »  0  (3-20) 

Equation  (3-20),  the  term  of  ASin  0t   is  not  always  equal  to  zero,  the  term 

of  SinG  should  be  zero.   Hence  the  value  of  8  will  be  either  zero  or  n. 

Inserting  equation  (3-7)  into  equation  (3-19): 

Ka  E  -.  3  Ka  E  3 
,  ,   ,       4    J  L .  1/2  (3-21) 

Ei-  F 

For  0=0.   and:  ~       g 

(KalEl,7  Ka3El  ,,*  (3.22) 

^  =        Et  +  F 

for  <9  *  n. 

As  a  check,  in  the  case  of  F  =  0,  i.e.  no  forcing  function  input,  equa- 
tions (3-21)  and  (3-22)  give  the  free  oscillation  frequency  of  system: 


10  = 


y'K  (at  +  f       a3E  )**  (3-23) 


For  the  linear  case,  i.e.,   a  =  0,  the  natural  frequency     : 

0J*   -  (axK)^  (3-24) 

It  is  exactly  the  same  as  we  discussed  in  the  linear  system. 

Recall  equations  (3-21)  and  (3-22);  the  value  of  lO     should  be  real, 
that  means: 

ai  +T  a3Ei2>  °  (3"25) 

for  both  the  phase  angle  is  either  zero  or  n.  And  : 

El    >  F  (3-26) 

for  the  phase  angle  is  zero  only. 

Equations  (3-25)  and  (3-26)  are  the  conditions  for  existence  of  the 
fundamental  frequency  forced  oscillation. 

A  free  oscillation  frequency  response,  i.e.   F  s  0  can  be  sketched 
from  equation  (3-23),  any  value  of  force  function  other  than  zero  also  can 
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be  sketched  on  both  sides  of  the  response  curves  of  F  ■  0.   It  Is  to  be 
noted,  that  the  phase  angle  between  the  force  function  and  the  error  signal 
is  opposite,  i.e.,  9  =  n,  when  the  response  curves  are  to  the  left  of  curve 
for  F  =  0,  and  the  phase  angle  between  them  is  in  phase,  when  the  response 
curves  are  to  the  right  of  it.   In  other  words,  it  is  according  to  whether 
the  frequency  is  less  or  greater  than  the  frequency  of  free  oscillation  for 
that  particular  amplitude  of  error  signal  with  a  constant  of  force  function 
input. 

A  typical  error  vs  frequency  (E. aJ  cu   )  response  curves  and  phase  re- 
lation sketch  for  a  different  value  of  forcing  function  and  different  char- 
acteristic of  nonlinearity  are  indicated  schematically  in  Fig.  3-5.   The 
response  curve  for  free  oscillation,  corresponding  F=0,  is  drawn  as  a 
dash  line. 

In  this  respect  the  behavior  of  the  nonlinear  oscillation  is  the  same  as 
that  of  the  linear  oscillation.   One  sees  that  the  response  curve  in  the 
nonlinear  cases  could  be  thought  of  as  arising  from  those  for  the  linear  case 
by  bending  the  latter  to  the  right  for  a  hard  spring  saturation  nonlinear- 
ity, and  to  the  left  for  a  soft  spring  saturation  nonlinearity. 

A  response  curve  in  "Input  vs  Frequency"  plane  (9  "^      cu       )  can  be 
also  obtained  from  relation  of: 

9c(t)   =  9r(t)  -  e(t)  (3-27) 

or: 

9  (t)   «  F  Cos(  uJ   t  +  9)  -  E  Cosh  t  (3-28) 

In  the  case  of  9  =  0,  equation  (3-28)  becomes: 


Hence 


©c(t)  =  (F  -  E  )  Cos^t  (3-29) 


9c  «  F  -  Ex  (3-30) 
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or; 

E     -  F  -   8  (3-31) 

Inserting  equation   (3-31)    into  equation   (3-21): 

.«*!    (6Q-F)     +  1     fa3   (6c-F>3,l/2 

02  *      (  q — '  (3-32) 

c 

For  the  case  of  8  «  if,  equation  (3-28)  becomes: 

8  (t)  *  (F  +  E.)  Cos'^t  (3-33) 

c  l 


Hence: 


Ej  »  -  (F  +  8c)  (3-34) 


Inserting  equation  (3-34)  into  equation  (3-23): 

K  «t(9  *F)  +  j  a3  (8,  +  F)3 

■  i 

C 

From  equations  (3-32  and  (3-35),  if  F  =  0,  the  result  is  the  same  form 
as  equation  (3-23),  that  means  the  response  curve  for  F  *  0  is  the  same  as 
plotted  in  the  E.  ■>-  Vd   plane,  only  the  difference  is  that  the  value  of  E.. 
changes  to  8  .  The  response  curves  for  the  values  other  than  F  ■  0  are 
changed,  and  the  phase  relationship  between  the  input  force  function  and  the 
output  is  reversed. 

There  is  a  numerical  example  for  plotting  the  response  curve  in  E.-^d) 
and  8  *v  oO    plane  respectively  in  section  3-7  of  this  chapter,  and  also 
experimental  response  curves  are  plotted. 

From  the  analysis  of  the  above  and  the  numerical  example  in  section  3-7, 
it  is  shown  that  a  saturation  nonlinearity  feedback  control  system  can  have 
fundamental  frequency  oscillations.  This  not  only  depends  on  the  character- 
istic of  the  nonlinearity,  but  also  depends  on  the  amplitude  of  input  force 
function. 
3-4  Forced  Oscillations:   Saturation  Nonlinearity  with  Damping. 

As  explained  previously,  if  there  is  no  damping  present  in  the  system, 
there  is  an  oscillation  error  signal  e(t)  ■  E.Cos(o)t)  either  in  phase  or 
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out  of  phase  of  180  degrees  with  the  forcing  function  of  Input.  In  the 
case  where  damping  Is  present  however,  the  forced  oscillation  displacement 
and  the  forcing  input  can  be  expected  to  be  out  of  phase,  just  as  In  the 
corresponding  of  linear  system. 

Recall  the  basic  equation  (3-6)  and  the  restoring  force  function  equa- 
tion (2-28): 

E  +/E  +  Kf(E)  «  -A  Cos('<  t  +  6)  -  BSin(  t  +  0)     (3-6) 
f(E)  =  axE  +  a3E3  (2-28) 

Then  equation  (3-6)  becomes: 

E  +  ^E  +  Ka  E  +  KaE3  «  -ACos(vt  +  G>  -  BSin(  t  -  0)    (3-36) 
Recall  one  solution  of  equation  (3-36)  is: 

e(t)  =  E.  Cos^t  (3-15) 

Inserting  equation  (3-15)  into  equation  (3-36): 

(-*.<.-    E,    +  Ka^,    +  7  Ka0E13)Cosr<  t   -o£u  E.Sin u  t  +  7  Ka.E.   Cos3     t 
1  11431  1  431 

-  (-ACos©  -  BSinG)  Cos  t  +  (ASin0  -B  Cose)  3in  t     (3-37) 
Neglect  the  higher  order  terms  of  harmonic  and  equates  the  coefficient 
of  Sin  u   t  and  Cos  t  respectively,  hence: 

(-^2E1  +  Ka1E1  +  I  Ka3Ex3)  =  ACosG  -  BSinS  (3-38) 

-«VE  «  A  Sine  -  B  Cose  (3-39) 

Squaring  equations  (3-38)  and  (3-39),  and  adding: 
(-  \:"El   +  Ka^  +  -  Ka3E13)2  +oc*tf\*   -  A2  +  B2  (3-40) 


or: 


(  EL  +  Ka1E1  +|  Kb  E  3)2  +*£    0)       Ex     s        'f2-      '    F  (3-41) 


also: 


22  .  <■       7  ^  2  2         2 

(El     -   F  )    -    fl,"L2Ej    (Ra1  +  f  K^   )    ~-  ~(^x  F  ) 

3 
+  Ej2    (Kat  +  4     KaE2)2     ^  0  (3-42) 
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As  a  check  of  equations  (3-38),  (3-39)  and  (3-41);  first  consider 
the  response  curve  in  the  high  frequency  and  the  low  frequency  ranges, 
the  phase  relation  between  the  error  signal  or  control  signal  and  the  in- 
put forcing  function  is  either  zero  or  180  degrees.   In  this  case,  the 
sine  terms  of  equation  (3-38)  and  (3-39)  are  zero,  the  response  curve  is 
the  same  as  equation  (3-20).   That  means  for  a  damped  system  if  the  phase 
angle  between  the  input  and  the  output  is  either  zero  or  180  degrees,  the 
frequency  response  is  exactly  the  same  as  the  system  with  no  damping. 

On  the  other  hand,  if  the  value  of  damping  is  very  small,  the  re- 
sponse curve  will  be  very  close  to  the  sketch  of  those  of  Fig.  3-5.   The 
only  difference  when  very  small  damping  is  present  is  that  the  response 
curves  are  rounded  off  in  the  vicinity  of  curve  for  F  *  0.  A  sketch  of 
them  is  shown  in  Fig.  3-6. 

For  the  phase  investigation;  recall  equation  (3-39)  and  rewrite  in 
this  form: 

SinO  -  /-  CosO  e^-E./F  (3-43) 

or: 

Sin(6  -  Tan"1 -    )  = 


F(   +   -)v"  (3-44) 

Hence,  the  phase  angle  between  the  error  signal  and  the  forcing  function 
of  input: 

'■'"''>(>,■■♦.■)   1/2     ♦Tan'1—  (3-45) 

Examine  equation  (3-45)  in  the  case  of  <-.'  ■  0,  6  is  either  zero  or  n, 
when  damping  is  present  the  value  of  6  is  a  function  of  :..  and  E-,  if  the 
values  of u   and  F  are  fixed.   It  should  be  noted,  there  are  two  values  of 
phase  shift  for  one  value  of  error  signal;  one  is  for  the  low  frequency 
(a)<.U}y\)   and  the  other  is  for  the  high  frequency  (cO^C\Jn);   one  is  from 
equation  (3-45)  and  the  other  i&  from  the  180  degrees  minus  the  value  of 
equation  (3-45). 
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For  the  existence  conditions  for  forced  oscillation  when  damping  is 
present,  recall  the  equation  (3-42),  and  put: 

Kax  +  -|  Ka3E1  -  X  (3-46) 

and: 

(Ex2  -  F2)     -(ZE^X  -v  (Ej2  -  F2)      +  EjV  -  0      (3-47) 
Solve  equation  (3-47): 


u)     '    ■  -  *<£&*)]*  J l(*e*x  -cf(e,*-Fzl  (3-48) 

It  is  to  be  noted  that  the  conditions  for  forced  oscillations  in  an 
undamped  system  are  still  valid  for  a  damped  system  within  the  low  and 
high  frequencies,  in  other  words,  it  is  valid  for  the  condition  of  phase 
angle  either  zero  or  «. 

From  equation  (3-48)  where  forced  oscillation  exists^  the  value  of 
should  be  real,  hence;  the  condition  for  equation  (3-48)  to  be  real: 


and; 


or: 


E2  -  F2>  0  (3-49) 


2E2X  -   '^(Ej2  -  F2)  .   0  (3-50 


x     __,'  (E  2  .  F2)   0  (3-51) 


[2Ej2X  -     *■*  (Ex2      -      F2)       j   2      -   4E12X2(El2   -   F2)   >  0  (3-52) 


222  22  ")    0  2? 

4EX   X  F       +     at*  (Ej      -   F   )      >     4Ej   X  (E^      -   F   )  (3-53) 


Equation  of  (3-49);  (3-51)  and  (3-53)  are  the  conditions  for  the 
forced  oscillation  of  a  damped  system. 
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and; 


or: 


3-5  Jump  Resonance: 

Recall  a  sketch  of  frequency  response  curve  of  a  soft  saturation  non- 
llnearlty  shown  in  Fig.  3-7.   When  a  experimental  frequency  response  is 
being  measured  for  a  system,  this  measurement  can  be  done  by  an  analog 
computer  to  measure  the  amplitude  of  the  input  and  the  output,  at  the  same 
time  measure  the  phase  relation  between  them. 


(9 


CO 


Fig.  3-7  The  Region  of  Jump  Resonance 

It  is  found  that  under  certain  conditions  of  amplitude  and  frequency 
an  infinitesmal  change  of  either  frequency  or  amplitude  of  the  input  signal 
causes  a  large  and  discontinuous  jump  in  amplitude  of  the  output,  at  the 
same  time,  a  discontinuous  phase  jump  occurs.   Fig.  3-8  shows  a  set  of 
amplitude  jump  and  phase  jump  of  a  typical  system  from  an  Electronic  Analog 
Computer  experimental  work. 

Fig.  3-7  is  shown  as  a  typical  system  with  a  soft  saturation  non- 
linearity  and  with  a  fixed  amplitude  of  the  input,  which  is  larger  than  the 
value  of  saturation  of  nonlinearity.   If  we  are  going  to  measure  the  frequency 
response,  being  with  a  low  frequency,  then  increase  the  frequency,  given  a 
response  from  A  to  B.  At  point  B  with  a  frequency  U),  ,   an  infinitesimal 
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increase  In  frequency  causes  the  amplitude  to  jump  upwards  from  point  B  to  C, 
and  with  a  phase  shift  from  below  90  degrees  to  a  value  above  90  degrees. 
Further  increases  in  the  frequency  cause  the  output  to  follow  the  gradual 
smooth  curve  to  D.   Now,  if  the  input  frequency  is  reduced  from  higher 
values,  the  output  will  follow  the  smooth  curve  DC.   When  the  frequency  is 
reduced  to  the  point  E,  if  the  frequency  is  infinitesimal  reduced,  the  out- 
put jumps  down  from  point  E  to  point  F. 

It  is  found  that  either  from  the  theory  analysis  or  from  the  experi- 
mental work,  a  hard  spring  saturation  nonlinearity  can  also  produce  the 
jump  resonance,  only  the  difference  between  them  is  the  jump  in  the  oppo- 
site direction  from  each  other.   A  typical  system  for  a  hard  spring  satura- 
tion nonlinearity  will  be  shown  in  the  experimental  section. 

From  the  above  discussion,  the  portion  of  the  curve  from  point  B  to 
point  E  cannot  be  obtained  experimentally,  usually  this  portion  of  the 
curve  corresponds  to  an  unstable  operating  condition. 
3-6  Lower  and  Upper  Jump  Frequency  of  the  Jump  Resonance: 

As  previously  discussed,  the  jump  frequency  of  (V,    is  called  "Upper 
Jump  Frequency"  and  the  jump  frequency  of  tU.^    is  called  "Lower  Jump  Fre- 
quency".  Recall  the  equation  (3-41)  and  the  sketch  of  the  response  curve 
shown  in  Fig.  3-7: 

(-d)  EL  +  Ka^j  +  ~     Ka3E13)2   +  <^V  E^   «<*/f  -+4A   F2     (3-41) 

The  equation  of  the  locus  of  the  vertical  tangents  can  be  found  by 
differentiating  equation  (3-41)  implicitly  with  respect  to  E.  and  setting 
dec     /dE.  equal  to  zero,  the  result  is: 

(-0)'     +  Kaj  +  |-Ka3Ex2)(    +  Ka][  +  -^-Ka^2)^       (3-54) 

Check  the  equation  of  (3-54),  if  the  damping  is  very  small,  the  term 


oC'oC       can  be  neglected,  that  leads  to  a  pair  of  equations 

ll  +  I  a3El 
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K(a.  +  ~  a  E.3)  =  0  (3-55) 


K(ax  +  -|  a3Ei2)  "  °  (3-56) 

The  curves  corresponding  to  these  equations  are  shown  in  Fig.  3-9. 
It  is  to  be  noted,  the  equation  of  (3-55)  is  the  same  as  the  response 
curve  equation  (3-23)  for  the  free  oscillation  of  an  undamped  system.   The 
equation  of  (3-56)  is  the  locus  of  the  vertical  tangents  of  curve  of  Fig. 
3-5. 

In  the  case  when  damping  is  present,  the  locus  of  vertical  tangents 
equation  (3-54)  will  be  directed  by  the  equations  of  (3-55)  and  (3-56), 
in  particular  there  should  be  one  branch  near  the  response  curve  for  free 
oscillation,  and  another  near  the  points  where  the  curves  for  the  undamped 
system.   In  other  words,  the  curve  must  appear  as  they  are  shown  in  Fig. 
3-10. 

Rearrange  the  equation  (3-54),  and  solve  for  tC     ,  which  leads; 

U3  ~\ii-ZK*,  r3Ka^;~W-)±  J(**W;  f   M  ;»;.:  (3-57) 


or: 


In  which      =  Ka.,  it  is  the  natural  frequency  of  the  linear  system.   If 
we  know  the  constants  of  system,  and  the  amplitude  of  fundamental  fre- 
quency forced  oscillation,  the  upper  and  the  lower  jump  frequency  can  be 
calculated  from  either  one  of  the  equations  (3-54),  (3-57),  and  (3-58). 


3-7  Analog  Computer  Analysis  for  a  Second  Order  System  with  Soft  Saturation 
Nonlinearity : 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  3-11. 
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Fig.  3-11  Block  Diagram  of  a  2nd  Order  System 
with  a  Soft  Saturation  Nonlinearity 

In  which  the  restoring  force  function  of  nonlinearity  is: 

f(E)  =2E-  0.032E3  (3-59) 

The  characteristic  curve  of  restoring  force  function  is  drawn  in  Fig.  3-12 

where : 

E  -  1.00   1.50  2.00  2.50  3.00  3.50  4.00  5,00  6.00 

f(E)«  1.97  2.89  3.75  4.50  5.15  5.65  5.95  6.00  5o10 

Before  the  experimental  analysis  consider  the  theoretical  calculations 

and  curves.   First  consider  the  response  curve  on  the  &,      .plane  for  the 

system  without  damping,  and  with  an  input  function: 

0  (t)  «  3  Cos(iit  +  6)  (3-60) 


Recall  the  equations  of  (3-21);  (3-22)  and  (3-23)  and  leads 

2  1/2 
a)   -  (4  -  0.048E  )  ' 


(3-61) 


for  free  oscillation  and: 


I 


4Et  -  0.048EJ 


1/2 


E  +  3 


(3-62) 


for  the  phase  angle  9  ■  n,  and: 


-  ( 
for  the  phase  angle  0=0. 


4E   -  0.048E.   1/2 


E  -  3 


38 


(3-63) 


ftE) 


-b.c 


6-.0 


rI^  3---  Sb  -r  eteristic  Curve  of  Restoring 
Function;  f(E)  =*  2E  -  0.032^ 
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The  calculations  and  curves  of  equation  (3-61)  to  (3-63)  are  shown 
in  Table  3-1  and  Fig.  3-13  respectively. 

Table  3-1:  Calculation  for  response  curves  on  the  E.^'  '■'  plane: 


Amplitude  of 

Frequency  for 

Frequency  for 

Frequency  for 

Error  Signal 

Free  Oscill. 

©  =  n 

9-0 

(   ) 

(   ) 

(   ) 

(   ) 

0.00 

2.00 

/ 

/ 

0.50 

1.996 

0.75 

/ 

1.00 

1.99 

0.99 

/ 

2.00 

1.95 

1.23 

/ 

3.00 

1.89 

1.37 

/ 

3.50 

l.K 

1.35 

4.88 

4.00 

1.80 

1.36 

3.60 

5.00 

1.68 

1.32 

2.65 

6.00 

1.53 

1.23 

2.16 

7.00 

1.29 

1.08 

1.73 

8.00 

1.00 

0.80 

1.22 

9.00 

0.35 

0.28 

0.42 

For  the  calculation  of  response  curves  and  phase  shift,  when  damping  is 
present,  recall  equation  (3-42),  and  leads: 

ZEl  E  2 

r  ).  0.252-   — 2 — p  <4  "  0.048EJ2)]   +     2*  2  (4  -0.048E2)2  =  0 
E  -3  E   -  3 

(3-64) 

for  the  phase  shift,  recall  equation  (3-45): 

e-Sln_1   F(^.^)1/2   +T«"'1-g-  <3-45> 

It  is  to  be  noted,  from  equation  (3-45),  and  the  response  curve  for  an  un- 
damped system  shown  on  the  e,  -vftj  plane,  when  the  frequency  is  larger  than 
the  frequency  of  free  oscillation  with  a  particular  amplitude  of  input, 
the  phase  angle  ©  is  zero  and  is  tt,  when  the  frequency  is  below  the  free 
oscillation  frequency.   Hence;  the  phase  calculation  from  equation  (3-45) 

should  be  modified  as: 

.   °'25E1  .  0,15 

0  »  180  -  (Sin   s-  n   +  Tan  (3=65) 

3(  *>\     0.25  y' 
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o 


o 


o 

V5 


Ui 


-lL 


for         ,  and: 

0.25E,         T   -1^ 

e  =  sm"1  4—  +  fan  (3-66) 

3(   +  0.25  y 
for 

If  the  response  curve  Is  presented  on  the      plane,  the  condition 
of  equation  (3-65)  and  (3-66)  should  be  reversed. 

The  calculation  for  the  system  with  damping  of  0.25  and  the  curves  are 
shown  in  Table  3-2  and  Fig.  3-14  respectively. 

The  jump  frequencies  of  system  can  be  calculated  from  equation  (3-57) 
and  also  can  be  seen  from  the  response  curve.   It  is  shown  schematically 
in  the  diagram  of  Fig.  3-14. 

Table  3-2:   Calculation  for  Response  Curves  (with  a  damping  of  0.25) 


Amplitude  of 

Frequency 

Phase   for 

Frequency 

Phase   for 

Error  Signal 

for  ui  <ti)n 

v)   s  tVri 

for  - 

•tX)^ 

(  &,  ) 

(  «   ) 

(  ?  ) 

(    CO) 

(        ) 

1.00 

1.00 

164.00 

1 

/ 

2.00 

1.24 

161.00 

1 

/ 

3.00 

1.35 

158.50 

1 

/ 

3.50 

1.36 

156.50 

4.83 

3.00 

4.00 

1.37 

155.00 

3.58 

7.00 

4.50 

1.40 

153.00 

2.80 

12.00 

5.00 

1.36 

151.50 

2.50 

15.50 

6.00 

1.24 

144.50 

2.05 

21.00 

7.00 

1.09 

135.00 

1.64 

28.50 

8.00 

0.87 

121.50 

1.16 

45.50 

8.70 

0.65 

90.00 

0.65 

90.00 

For  the  experimental  analysis,  the  first  investigating  is  the  restor- 
ing forcing  function.   It  is  to  be  approximated  as  a  two  slope  lines,  one 
slope  is  2,  and  the  other  is  zero.   The  circuit  setting  in  the  analog  com- 
puter is  shown  below,  and  also  all  components  to  be  used  are  schematically 
in  the  circuit  of  Fig.  3-15. 

The  characteristic  curve  for  the  restoring  forcing  function  from  the 
computer  is  shown  in  Fig.  3-16;  (a)  is  shown  the  waveform  of  input  (e  ) 
and  output  (e  ),  (b)  is  shown  the  characteristic  of  the  output  to  the  input. 
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Fig.  3-15  Operational  Diagram  of  a  Function 


Generator   for     f(E) 


2E 


0.032E' 


The  operational  diagram  of  analog  computer  setting  up  for  system  of 
block  diagram  in  Fig.  3-11  is  shown  in  Fig.  3-17,  in  which: 


e(s)  = 


<(? 


oce. 


c(s)  =  I^± L 


l(s)  -  0t  (s)J 


(s) 
(s) 


(3-67) 
(3-68) 
(3-69) 


From  equations  (3-67)  to  (3-69)  system  scaling  is  shown  below: 
where:     e^o  m   1.00;  tU.m   1.00;  ^c«  1.00;  *&  m   2.00 


Wl  *  al  Rfl/Rl  =  1,0° 
W2  =  a2  Rfl/R2  =  1,0° 
W3  *  Rf2Cfl/a,  =  4.00 


1 


1.00 
1.00 

0.50 


W4  "  a3Rf2/a4R3  *  4*00 

a3  =   1.00; 

R3  *  0.50; 

W5   =  a5/   R5Cf2   «  2.00 

a5  =   1.00; 

R     »  0.50 

cf2  =  1.00 

W6  "  a6Rf     /R6  =   U0° 

a,   =   1.00, 
o 

Rft  -   1.00 

Rf3  -   1.00 

lfl 


1.00; 


R2  =  1.00 


vf2 


1.00 


in  which  R  in  megohm;  C  in  U  f 


RL  m   1.00 
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The  results  ftom  the  analog  computer  experimental  work  Is  condensed 
in  Table  3-3;  and  the  response  curve  with  phase  shift  represented  on  i, 
plane  is  shown  in  Fig.  3-18  respectively.   Also  the  jump  frequencies  are 
indicated  schematically  in  these  diagrams. 

A  theoretically  calculated  response  curve  in  /:,  ,L:   plane  for  comparing 
is  shown  with  the  experimental  response  curve  of  Fig.  3-18.   Comparing  the 
results  from  the  computer  experimental  work  with  the  results  from  theoretical 
calculation,  it  is  shown  that  the  results  are  very  close  within  high  fre- 
quency and  low  frequency  range.  There  is  a  resonance  phenomena  between  the 
jump  frequencies. 
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Table   3-3  Data    for  Response  Curves    from  Analog  Computer   for  System 

»  0.25;   K  «  2;    f(E)   =  2E  -   0.032E3;        ,'.'.  (t)   «  3.0Cos(,t   +  0) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Tnput 

Error  Signal 

Error 

Output 

Output 

((X)   ) 

(   ) 

(  -  ) 

(   ) 

(  -  ) 

0.500 

0.38 

180.00 

3.20 

0.00 

0.628 

0.50 

180.00 

3.40 

0.00 

0.755 

0.75 

175.00 

3.60 

5.00 

0.880 

1.25 

168.00 

4.00 

9.00 

1.000 

1.50 

163.00 

4.50 

12.00 

1.130 

2.15 

158.50 

5.00 

24.00 

1.260 

3.00 

150.00 

5.40 

35.00 

1 . 260<+) 

9.00 

/ 

6.50 

/ 

1.380 

8.25 

25.00 

5.50 

140.00 

1.510 

7.30 

14.00 

4.50 

155.00 

1.640 

7.00 

5.00 

3.80 

159.00 

1.760 

6.50 

0.00 

3.40 

165.00 

1.890 

6.20 

0.00 

3.00 

172.00 

2.010 

6.00 

0.00 

2.60 

180.00 

2.140 

5.50 

0.00 

2.20 

180.00 

2.260 

5.20 

0.00 

2.00 

180.00 

2.390 

5.00 

0.00 

1.60 

180.00 

2.510 

4.80 

0.00 

1.40 

180.00 

2.830 

4.40 

0.00 

1.20 

180.00 

3.140 

4.00 

0.00 

1.00 

180.00 

3.460 

3,80 

0.00 

0.80 

180.00 

3.770 

3.70 

0.00 

0.60 

180.00 

4.090 

3.60 

0.00 

0.50 

180.00 

(b)  Frequency  Decreasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(.0  ) 

(  B,) 

(   ) 

(  c   ) 

(  *   ) 

2.510 

4.75 

0.00 

1.50 

180.00 

2.200 

5.20 

0.00 

2.00 

180.00 

1.890 

6.20 

0.00 

3.00 

180.00 

1.760 

6.80 

0.00 

3.20 

173.00 

1.640 

7.10 

0.00 

3.80 

168.00 

1.510 

7.58 

o.OO 

4.50 

162.00 

1.380 

8.30 

12.00 

5.50 

156.00 

1.260 

9.20 

23.00 

6.80 

140.00 

1.130 

10.00 

34.00 

8.00 

135*00 

1.000 

11.50 

65.00 

9.50 

120.00 

0.880 

11.00 

/ 

8.70 

/ 

0.880(-) 

0.75 

175.00 

4.00 

0.00 
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The  resonance  phenomena  in  Fig.  3-18  is  actually  similar  to  the  re- 
sonance in  the  linear  system.   It  depends  on  the  value  of  damping  in 
the  system.   Another  experimental  result  for  the  same  system  and  same 
value  of  input,  only  changed  the  damping  from  0.25  to  0.5,  and  is  shown 
in  Table  3-4  and  Fig.  3-19  respectively.   In  this  case  the  resonance  pheno- 
mena is  much  lower  than  the  case  of  a  damping  0.25. 

It  is  to  be  noted,  if  the  damping  of  a  system  is  increased  the  range 
of  jump  frequencies  is  decreased.   For  this  reason,  when  the  damping  in- 
creased to  limit  value,  at  which  no  jump  phenomena  existed,  it  means  that 
the  system  becomes  a  linear  system.   On  the  other  hand,  if  the  system  is  an 
undamped  or  a  very  low  damping  system,  it  will  be  very  oscillatory  or  a 
limit  cycle  exists. 

The  response  curve  represented  on  "Output  vs  Frequency"  (<P<.-\  c;  ) 

plane  can  be  calculated  from  equations  (3-32)  and  (3-35).   An  experimental 

response  curve  shown  on  the  6>c  >  (fi   plane  from  the  system  with  a  damping  of 

0.25  is  shown  in  Fig.  3-20. 

3-8  Analog  Computer  Analysis  for  a  Second  Order  System  with  Hard  Saturation 
Non linearity: 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  3-21:  in  which 


0rO) 


v^g) 


•:> 


— \~> 


/.-s 


SCS  1 


Qc(S) 


Fig.  3-21  Block  diagram  of  a  Second  Order  System 
with  a  Hard  Saturation  Nonlinearity 

assume  the  restoring  forcing  function  of  the  nonlinearity  is: 

f(E)  =  0.425E  +  0.005E3 


(3-70) 
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Table  3-4,   Data  for  Response  Curve  from  Analog  Computer  for  System 

U    «  0.50;  K  «  2;   f(E)  -  2E-0.032E3;   6,  (t)  *  3.0Cos(  < ' t+e) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(-  ) 

(  e,) 

(  -  ) 

(  o 

(  6  ) 

0.628 

0.50 

173.00 

3.40 

7.30 

0.755 

0.75 

165.00 

3.50 

12.00 

0.880 

1.00 

155.00 

4.00 

16.50 

1.000 

1.50 

140.00 

4.75 

24.00 

1.130 

1.75 

125.00 

6.00 

30.00 

1.260 

8.30 

91.00 

8.00 

85.00 

1.380 

8.10 

47.00 

6.25 

120.00 

1.510 

7.50 

36.00 

5.25 

135.00 

1.640 

7.00 

25.00 

4.25 

148.00 

1.760 

6.50 

20.00 

3.50 

160.00 

1.890 

6.00 

12.00 

2.75 

172.00 

2.010 

5.50 

5.00 

/ 

/ 

2.140 

5.00 

0.00 

/ 

/ 

2.260 

4.50 

0.00  . 

1.75 

180.00 

2.390 

4.25 

0.00 

/ 

/ 

2.510 

4.00 

0.00 

1.25 

180.00 

2.810 

3.75 

0.00 

1.00 

180.00 

3.140 

3.50 

0.00 

0.75 

180.00 

3.460 

3.35 

0.00 

0.50 

180.00 

3.770 

3.25 

0.00 

0.30 

180.00 

4.090 

3.10 

0.00 

/ 

/ 

(b)  Frequency  Decreasing 


equency  ol 

:   Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(  a)  ) 

(  /  ) 

(  *  ) 

(  <S>t) 

(  L     ) 

2.510 

4.30 

0.00 

1.25 

180.00 

2.200 

5.00 

0.00 

1.75 

180.00 

1.890 

5.75 

0.00 

2.75 

172.00 

1.760 

6.25 

7.50 

3.50 

160.00 

1.640 

7.00 

15.00 

4.10 

148.00 

1.510 

7.75 

27.00 

5.00 

132.00 

1.380 

8.25 

32.00 

6.50 

120.00 

1.260 

8.50 

60.00 

7.50 

97.00 

1.130 

7.00 

80.00 

8.20 

40.00 

1.000 

2.00 

130.00 

4.75 

19.50 

0.880 

1.25 

145.00 

4.00 

0.00 
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Fig,   3-23£Charact<3ristic  Curve  oi" 
Restoring  faction:  f(^)-0.425r*0.005E3 
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The  characteristic  curve  of  equation  (3-70)  Is  shown  in  Fig.  3-2la,  where: 
E  -  1.00;   2.00;   3.00;   4.00;   5.00;   6.00;   7.00;   8.00. 
f(E)  =  0.43;   0.89;   1.42;   2.02;   2.75;   3.63;   4.69;   5.96. 

A  theoretical  analysis  can  be  calculated  from  equations  (3-70)  and 
(3-21)  to  (3-23),  as  we  did  in  section  3-7.   Only  the  analog  computer  anal- 
ysis is  considered  in  this  section. 

For  the  simulation  of  the  restoring  forcing  function  of  nonlinearity 

\ 
from  the  analog  computer,  it  is  approximated  by  two  straight  lines,  one 

with  a  slope  of  0.425,  and  the  other  with  a  slope  of  2.00.   The  circuit  set- 
ting in  the  computer  is  shown  in  Fig.  3-22,  and  all  components  to  be  used 
are  schematically  in  the  diagram. 


<£.*■■* 


A  C 


- 


-  J>i) 


Fig.  3-22  Operational  diagram  of  the  Function  Generator 
for   f(E)  =  0.425E  +  0.005E3 

The  characteristic  curve  of  restoring  forcing  function  from  the  com- 
puter is  shown  in  Fig.  3-23. 

The  operational  diagram  of  analog  computer  set  up  for  the  whole 
system  is  mostly  the  same  with  the  diagram  of  Fig.  3-17,  the  only  dif- 
ference is  the  function  generator  circuit. 
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The  results  from  the  analog  computer  are  condensed  in  Table  3-5, 
and  the  response  curves  with  phase  shift  shown  on  the  £/•*  and 
plane  are  shown  in  Fig.  3-24  and  3-25  respectively.   In  which  the  input 
forcing  function  is  9      (t)  =3.50  Cos(<  t  +  9). 
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Table  3-5.   Data  for  Respor se  Curves  from  Computer  for  a  System  with  a 

Hard  Saturation  Nonl  inearity,  in  which:  oC     =  0.25;  K  »  2; 

f(E)  =  0.425E  +  0.005E3;  Q   (t)  =  3.5  Cos  (..it  +  G). 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(rt) 

(  e,) 

{0  ) 

(  ®c  ) 

(  (9  ) 

0.500 

1.75 

180.00 

5.50 

0.00 

0.628 

3.50 

180.00 

6.50 

4.00 

0.754 

6.50 

175.00 

9.00 

7.50 

0.880 

8.00 

168.00 

10.50 

10.00 

1.000 

9.00 

164.00 

12.00 

1.1.50 

1.130 

10.50 

158.50 

13.00 

13.00 

1.260 

12.25 

152.00 

14.50 

15.00 

1 .  380 

14.00 

145.50 

16.30 

17.00 

1.510 

16.25 

141.00 

18.75 

20.00 

1.640 

19.00 

134.00 

21.30 

24.00 

1.760 

22.25 

124.00 

22.00 

37.00 

1.890 

26.00 

102.00 

25.00 

45.00 

1.950 

26.50 

/ 

25.50 

55.00 

2.010 

26.00 

24.00 

25,00 

90.00 

2.140 

4.50 

20.00 

0.86 

150.00 

2.260 

4.35 

15.00 

0.65 

165.00 

2.510 

4.25 

7.50 

0.50 

180.00 

2.830 

4.15 

0.00 

0.40 

180.00 

3.140 

4.00 

0.00 

0. 10 

180  00 

^.770 

3.80 

0.00 

/ 

/ 

(b)  Frequency  Decreasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(CO  ) 

(  ^/) 

(  &  ) 

(        ) 

(   &  ) 

2.510 

4.25 

5.00 

0.50 

180.00 

2.400 

4.35 

iQtKr 

0.63 

180.00 

2.270 

4.50 

13.50 

0.68 

172.50 

2.140 

4.60 

15.00 

0.75 

170.00 

1.890 

4.85 

21.00 

1.13 

157.50 

1.760 

5.00 

23.50 

1.25 

155.00 

1.640 

5.35 

25.00 

1.60 

152.00 

1.510 

5.75 

26.00 

2.13 

150.00 

1.380 

6.75 

27.50 

16.30 

18.50 

1.320 

12.25 

/ 

/ 

/ 

1.260 

11.70 

150.00 

14.50 

15.50 

1.130 

10.25 

155.00 

13.00 

13.50 

1.000 

9,00 

160.00 

11.70 

12.00 

0.880 

7.50 

172.00 

10.50 

10.00 
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CHAPTER  IV 

SUB-HARMONIC  OSCILLATION  IN  NONLINEAR  FEEDBACK  CONTROL  SYSTEM 
4-1     General   Description? 

Up  to  the  preceding  chapter,  there  is  only  discussed  the  forced 
oscillations  for  which  the  frequency  is  the  same  as  that  of  the  external 
force.   A  nonlinear  control  system,  under  suitable  conditions  may  have  a 
steady  state  oscillation  in  which  the  main  component  has  a  frequency  which 
is  a  fraction  1/n  (n  is  a  real  and  positive  integer)  of  the  forcing  func- 
tion frequency.   These  oscillations  are  called  "Sub-Harmonic  Oscillation". 

Recall  a  linear  system  (Section  2-3),  there  are  two  oscillation  terms 
in  the  transient  response,  if  the  frequency  of  the  free  oscillation  is 
<*)  /n  (n  is  an  integer),  then  a  force  function  of  frequency  «j    can  excite 

the  free  oscillation  in  addition  to  t  he  forced  oscillation  of  frequency^  . 

-  t/2 
But  as  t  increases,  e       damps  the  free  oscillation  until  only  the 

steady  state  oscillation  remains.   Hence  there  is  no  sub-harmonic  or  super- 
harmonic  oscillation  existing  in  a  linear  control  system  at  steady  state. 

In  the  case  of  nonlinear  control  system,  the  steady  state  response 
may  not  have  the  same  frequency  of  oscillation  as  the  frequency  of  the 
forcing  function.   In  fact,  the  frequency  of  the  system  response  depends 
not  only  upon  the  frequency  of  forcing  function  but  also  upon  the  amplitude 
of  the  input.   The  principle  of  superposition  no  longer  holds  and  the  linear 
circuit  characteristic  is  not  applicable.   In  general,  it  can  be  said  that 
the  nonlinear  system  can  have  wide  variety  of  almost  periodic  oscillations, 
the  frequencies  of  which  differ  from  the  frequency  of  the  applied  force  and 
vary  with  time  as  well  as  with  different  initial  conditions. 

Sub-harmonic  oscillation  always  existed  with  a  frequency  higher  than 
the  frequency  of  free  oscillation  with  a  particular  amplitude  of  error  signal. 
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The  response  curve  which  is  a  part  of  an  ellipse  or  a  hyperbola  in  the 
E«~«  plane  is  always  nearly  parallel  to  the  right  side  of  the  free  oscil- 
lation response  curve.   A  typical  sketch  of  a  soft  saturation  nonlinearity 
is  shown  in  Fig.  4-1. 


&   >/j   <A       xA.         30b- liar  moHi'i. 


■ 


Fig.  4-1  Fundamental  and  Sub-harmonic  Response 

Curves  for  a  Soft  Saturation  Nonlinearity 

Under  some  suitable  conditions,  there  may  exist  any  order  of  sub- 
harmonic  oscillation,  or  only  one  type  of  order  can  exist  depending  on 
the  characteristic  of  the  nonlinearity  and  the  input  forcing  function. 

Sub  harmonic  oscillation  and  response  of  nonlinear  vibratory  systems 
are  often  discussed  in  the  field  of  mechanical  system  *  .  There  is  no 
body  of  knowledge  applied  to  feedback  control  servomechanism.   For  the  feed- 
back control  system,  there  is  a  standard  differential  equation  (3-6a),  it 
is  very  similar  to  the  equation  for  the  mechanical  system,  the  only  dif- 
ference is  that  the  forcing  function  of  a  feedback  control  system  is  a  func- 
tion of  input  frequency. 

Although  any  order  of  sub-harmonic  oscillation  can  exist  in  a  non- 
linear feedback  control  system,  in  this  paper  only  the  sub-harmonic  of  order 
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2  and  3  are  discussed. 

4-2  Sub-harmonic  oscillation  of  Order  2  without  Damping: 

Sub-harmonic  oscillation  of  order  2  can  exist  in  most  nonlinear  feed- 
back controls.   When  it  occurs,  the  frequency  of  output  (  C ,-. )  or  error 
(e  )   is  only  1/2  the  frequency  of  input  (  6Jr   ).   A  typical  diagram  from 
the  computer  experimental  work  shown  on  the  "Input  vs  Output"  plane  and  the 
waveforms  of  "Input  and  Output"  are  shown  in  Fig.  4-2  and  4-3,  respectively. 

Recall  the  basic  equation  for  forced  oscillation  in  nonlinear  system 

with  un-damping  Eq.  (3-11). 

E  +  Kf(E)  =  -A  Cos(  t£    t  +  ©)  (3-11) 

For  the  same  type  of  restoring  force  function  with  chapter  3,  leads: 

E  +  KajE  +  Ka3E3  =  -  A  Cos(^  t  +  6)  (3-12) 

Assume  one  solution  of  equation  (3-12)  is: 

E  =  E  2Cos  -jvOt   +  E1   Cos  to   t  (4-1) 

Where  E1  ,2   is  the  amplitude  of  sub-harmonic  of  order  2,  E.  is  the  amplitude 

of  fundamental  frequency,  since: 

-     1  1 

E  =  j   "  E..?  Cos  -jiCt   •«]  E.  Cos-  t  (4-2) 

(£&  +fl,       €u)t  f  2F£zH 

(4-3) 
We  are  only  interested  to  the  terms  of  second  sub-harmonic  and  fund- 
amental frequency,  inserting  equations  (4-1)  to  (4-3)  into  equation  (3-12) 
the  result  is: 

(-  $rs£.  +  W/%  +  £  K*3*jL+£  <%  £y,  E?  )  Cc'-  '-t  Wt   t  (-  ^>V;  - 

(4-4) 
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■:-3  2nd  Order  Sub-harmonic  waveforms 
for  c.  System:   f(B)  =2E-0,0048e3 
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As  dtscussed  In  Chapter  3,  Section  3-3,  if  the  frequency  is  larger 
then  the  frequency  of  free  oscillation,  the  phase  angle  0  is  zero  on  the 
"Error  vs  Frequency"  plane,  hence  equation  (4-4)  leads: 

-  4-t.  i    6  ^a3e,^  (4-5) 


C°     ~~  -F  (4-6) 

Equation  (4-5)  is  the  condition  for  existence  of  the  sub-harmonic  of  order 
2.   Solve  E  /2: 

Since  E^  .«  should  be  real;  which  leads: 


00 


(4-Ka,  +  /S/rajerJ^-  a3  %  o  (*-8) 


For  the  first  iteration  method,  consider  a  equal  to  zero: 
Equations  (4-5)  and  (4-6)  leads: 

O)-  ~<4K(>-<   ~4i  (4-9) 

t,   *  £f  (4-io) 

Inserting  equations    (4-9)   and   (4-10)    into  equation   (4-8): 

CO  ^   2>(a)a  t  8  ifa3  Fx)&  it3  <   0  (4-11) 

From  Equation  (4-5),  it  represents  an  ellipse  or  a  phperbola  in  the  E.»~ 

vs  u>  plane  depending  on  the  sign  of  a  a_,  also  w    has  a  maximum  when  a_ 

■i  j 

is  less  than  zero,  and  a  minimum  when  a,  is  greater  than  zero.  Actually, 

the  sub-harmonic  oscillation  only  existed  within  an  interval  of  frequency. 

That  means  only  one  part  of  ellipse  or  hyperbola  can  exist,  as  shown  in  Fig. 

(4-4). 

4-3   Sub-harmonic  Oscillation  of  Order  2  with  Damping: 

Recall  the  basic  equation  (3-41)  for  a  damped  system: 

E  +  <<^E  +  Ka.E  +  Ka3E3  =  -  ACos(^  t  +  9)  -  BSin(<t  t  +  0)        (3-41) 


66 


Assume  one   solution  of  equation    (3-41)    is: 

B —   :                             i    Elf,  (4- 
Since: 

U<A**OJt  r  £,            (Ob  (4-13) 

B          ;                                  £ifik             b  -Bi&$>»»&  (4-14) 

£  Mj&*        /''    '  '■■    ' 

r(^rBl          -4    T£fyBm)^»Mb  (4_15) 

Inserting  equations  (4-12)  to  (4-15)  into  equation  (3-41)  and  equal- 
ing the  coefficients  of  Cos  -xCUt   ,  Coso>£,  and  Sink'ir,  give  a  result  of 

•~^-%  n                      '                ;;  — o  (4-i6) 

=  -/*  <?<?$  £>  -  #  ^  #  (4-17) 

=  /}^  #  ~  £5  Cpi  (P  (4-18) 


From  equation  (4-18),  it  is  the  coefficient  of  the  term  of  Cos  -r  &  t, 


whe 


re     E..2  ^y hence* 
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(l)  Fundamental  Frequ 
( " )  2n  "  harmonic . 
(3)  5~cl  Sub  harm  nic. 


(a)  Soft  Saturation;  a^  -  0 


(b)  H^rd  i      m  <-~o       Q 


Fig  4-4  Response  Curves  for  ■ 

Sub  h:.r~crn.c  Oscillations 
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(4-19) 


Solve   for  E. .~: 


E 


1/2   '         '  (4-20) 


The  value  of  E. /2   should  be  real   therefore,    the  existence  of  2nd  order 
sub-harmonic  with  damping: 


.       .     -                                   a3   <     0  (4-21) 

From  equations  (4-16)  to  (4-18)  for  the  first  iteration  method,  for 
a_  equal  to  zero:  hence: 

*  4Ka  (4-22) 

(Ka   -     )E    +  ad       E   »  -    FCos©  -  .:-  FSin©  (4-23) 

(Kax  -  vJ   )E   -  *od  E.    ■  ^  FSin©  -  **   FCos©  (4-24) 
Solve  E]A  and  E.   from  equations  of  (4-22)  to  (4-24): 


E1A  "  fFCC'         ^~pr  '  ]  (4=25) 


E    _ 2jC*).&*.  (4-26) 

IB 


Inserting  equations  (4-25)  and  (4-26)  into  equation  (4-21).  A 
general  equation  for  existence  of  2nd  order  sub-harmonic  oscillation 
with  damping: 


+  0Ka't    I-  "J  J  (4-27) 


For  the  case  of  8  equal  to  0  or  n,  i.e.,  Sin©  equal  to  zero  equation 
(4-27)  becomes: 


*,% 


0)  (A-28) 
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From  equations  (4-11)  and  (4-28)  they  were  shown  that  either  the 
system  with  damping  or  without  it  the  second  order  suh-harmonic  oscilla- 
tion only  existed  with  a  frequency  less  than  the  frequency  of  two  times 
the  natural  frequency  of  linear  system  for  the  case  of  a_  less  than  zero, 
and  greater  than  this  frequency  for  the  case  of  a.  greater  than  zero. 
4-4   Sub-harmonic  Oscillations  of  Order  3  without  Damping: 

Typical  waveforms  with  jump  phenomena  of  3rd  order  sub-harmonic 
forced  oscillation  is  shown  in  Fig.  4-5;  (a)   shows  the  transient  from 
the  2nd  order  to  the  3rd  order  sub-harmonic;  (b)  is  shown  reversing. 

For  investigating  the  existence  conditions  of  3rd  sub-harmonic  forc- 
ed oscillation  recall  the  basic  equation  of  (3-11);  and  leads: 

E   +  KajE  +  Ka  E3  «  -  A  Cos(<*  t  +  G)  (3-12) 

Assume  one  solution  of  equation  (3-12)  is: 

E  *  E./3  Cos  -  A    t  4  E.  Cos/v  t  (4-29) 


in  which  E, ,,   is   the  amplitude  of   3rd  order  sub-harmonic  oscillation 
hence: 


E  -  -     -~     E./3  Cos  -r    .    t   -  t>.   E     Cosivt  (4-30) 

E  "  *  E1/3(e21/3  +  E1/3E1  +  2E12'  COS3  *   * 

3 
|  (E1/33  +  6E1/3"  Ex  +  3E1  )  Cos   t  +.....  (4-31) 

From  equation  (4-29)  to  (4-31)  and  equation  (3-12)  leads. 


IK  (*-32) 

i'A  (4-33) 

Equation  (4-32)  E.  .«  is  not  equal  to  zero,  hence: 

(4-34) 
Solve  E...-  from  equation  (3-34). 
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Fig. h-S  Waveforms  of  Upper  and  Lower  Jimp  Transient 
for  .'•-  Typical  3rd  Order  Sub-harmonic 
"breed  Oscillation 
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(4-35) 


The  value  of  E..~  must  be  real,  then: 


-*7£>*  (4-36) 

Equation  (4-36)  is  the  condition  for  the  existenceof  sub-harmonic 
oscillation  of  order  ?   of  o  without   damping  system: 

From  equations  (4-32)  and  (4-33),  by  the  first  iteration  method, 


for  a  «  0;  then: 


-  9  Ke,  (4-37) 

(K&1   -    a     )Ej  *    F  (4-38) 

Solve  E   froir  equation  (4-37)  and  (4-38): 

Ej  -  |  F  (4-39) 

Inserting  equation  (4-39)  in  equation  (4-36)  the  condition  for  a  3rd 
sub-harmonic  oscillation  becomes: 


a  *l  a  v  i   -  (4-40) 

Also  the  response  equation  of  equation  4--3-1  becomes: 

ej     :  0  (4-4D 

The  equation  of  (4-41)  represents  an  ellipse  or  a  hyperbola  in 
the  E.  .»  vs  ''      plane  depending  on  the  sign  of  a  ,  also    has  a  maximum 
when  a-  0,  and  a  minimum  when   a_  >  0.  A  sketch  of  a  3rd  sub- 

harmonic  forced  oscillation  is  shown  in  Fig.  4=4. 
4-5  Sub-harmonic  Oscillation  of  Order  3  with  Damping: 
Recall  equation  (3-41): 

E   +  v'  E  +  KaxE  +  Ka3E3  =  -  ACos('t  +  G)   -  BSln(fc)  t  ♦  Q)      (3-41) 
Assume  one  solution  of  equation  (3-41)  is: 

E  -  E1/3Cos  -  «A)  t  +  E1ACos^t  +  EjgSin.  c  (4-42) 
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Since: 

E  «  E.  ._  •'"  Sin  \  •<   t      -  E1A  «t  Sin  .  t  +  Eln  «  Cos   t  (4 

1/3  -     3         IA  IB 

1 
E  ■  "  «Ei  /o  Cos  T  u  t  "  E1A  -  Cos   t  -  -  E.n   Sin  t  (4-44) 

9  1/J      3         IA  IB 

g3  "  !  El/3(El/32  +  El/3E1A  +  2E1A2  +  2Eib2)CoS  3^  + 

!E1/32E1BS14  *  t  +  1/4(El/33  *  6E1/32e1A  +  3E1A3  + 

3E1AE1B2)   Cos/-'  t   +  1   E1B  (E1B2  +  2  E1/32  +  i  2)  Sin*,  t  +  ..(4-45) 

4 

From  equation  (4-42)  to  (4-44)  and  (3-41),  by  equaling  the  coefficients 
of  Cos  —   t,  Sin  t,  Cos  t,  leads: 

(  K  I  (4-46) 

(4-47) 

(fat  (4-48) 

Rearrange  equation  (4-46)  and  solve  it  for  E.  .^   in  terms  of    , 

E1A»  and  E1B:  

By,    ■  (4-49) 

here  E1  »_  should  be  real,  hence: 

)]  (4-50) 

or: 

60   ^  .-  -,  (4-51) 

Equation  (4-51)  is  the  condition  of  3rd  sub-harmonic:  oscillation 

existing  in  saturation  nonlinearity  with  damping  system. 

By  the  iteration  method,  equation  (4-46)  to  (4-48)  leads: 

-  9Kax  (4-53) 

(Ka  -  CO   )EJA  +  .-.'   ETB  -  -   FCos  9  -  <     FSin©       (4-54) 

(Ka   -   )E1B  -  *        E1A  -  ,.     FSin  6  -     FCos©       (4-55) 
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Solve  equation  (4-53)  to  (4-55), 

(4-56) 


E1A   ' 

' 

E1B- 

. 

(4-57) 

Substituting  equations  (4-56)  (4-57)  into  equation  (4-52)  and  be- 
comes in  a  general  form: 


$rti~iizr< 


'  ±±*±hL  (4.58) 


In  the  case  of  G  equal  to  0  or  jt,  i.e.  Sin0  equal  to  zero  equation  (4-58) 
becomes : 


(4-59) 
In  conclusion,  the  condition  for  existence  of  sub-harmonic  oscilla- 
tion in  a  nonlinear  feedback  control  system  depends  on  the  characterics 
of  the  nonlinear  element  and  the  amplitude  of  input.   The  response  curve 
of  any  order  sub-harmonic  oscillation  with  a  saturation  nonlinearity  is 
either  an  ellipse  or  a  hyperbola  which  depends  on  the  characteristic  of 
nonlinearity.   Also  the  frequency  of  a  sub-harmonic  oscillation  existed 
may  be  either  a  maximum  when  a  soft  saturation  nonlinearity  is  presented 
or  a  minimum  when  a  hard  saturation  nonlinearity  is  presented.   The  limit 
value  of  minimum  or  maximum  is  almost  equal  to  the  number  of  sub-harmonic 
order  times  the  natural  frequency  of  linear  system. 

4-6  Analog  Computer  Analysis  for  a  2nd  Order  Sub-harmonic  Oscillation 
System: 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  4-6 t 
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Fig.  4-6  Block  Diagram  of  a  2nd  Sub-harmonic 
Oscillation  System. 

Assume  the  restoring  force  function  of  nonlinearity  is: 

f(E)  =  2E  -  0.0048E3  (4=59) 

The  characteristic  curve  for  equation  (4-59)  is  shown  in  Fig.  4-7,  in 
which: 

E  »  1.00;   2.50;   5.00;   7.50;   10.00;   12.50;   15,00. 
f(E)  =  1.995;   4.295;   9.400;  12.97;  15.20;   15.40;   14.80. 

For  the  simulation  of  restoring  force  function  of  nonlinearity  from 
computer,  it  is  approximated  from  two  straight  lines,  one  of  which  with 
a  slope  of  2;  and  the  other  with  a  slope  of  zero.   The  circuit  setup  in  the 
computer  is  identical  with  the  circuit  shown  in  Fig.  3-15.   The  only  dif- 
ference is  the  value  of  bias  voltage  to  be  used.   The  characteristic  curve 
of  restoring  force  function  from  the  computer  is  shown  in  Fig.  4-8. 

An  operational  diagram  setup  in  the  analog  computer  is  almost  the 
same  with  the  diagram  of  Fig.  3-2C,  only  one  difference  is  that  the  value 
of  constant  to  be  used. 

The  results  from  the  computer  are  condensed  in  Table  4-1  and  Fig.  4-9. 
Only  the  response  curve  is  of  interest,  hence,  only  one  response  curve  in 
the  error  vs  frequency  plane  is  presented.   It  is  shown  that  there  are  two 
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jumps,  one  is  the  fundamental  frequency  Jump,  and  the  other  is  the  second 
order  sub-harmonic  frequency  jump.  These  jump  phenomena s  are  shown  sche- 
matically in  the  diagram  of  Fig.  4-9. 

Table  4-1:  Data  for  a  2nd  Order  Sub-harmonic  Oscillation  Response 
Curve,  in  which:  f(E)  «■  2E  -  0.0048E3;  K  =  4,«/«  0.25; 
and  0  (t)  «  2.5  Cos(  ai   t  +  0) 


(a)  Frequency  Increasing: 

Frequency  of   Amplitude  of  Amplitude  of 

Input       1st  Harmonic  2nd  Sub-harmonic 

(«  )         (  C'i)  <*J  ) 


Remarks 


0*628 

0.22 

0.754 

0.33 

0.880 

0.45 

1.000 

0.50 

1.130 

0.63 

1.260 

1.00 

1.380 

1.15 

1.510 

1.50 

1.640 

2.00 

1.760 

2.50 

1.890 

3.00 

1.950 

3.25 

1.980* 

10.50 

*lst  harmonic 

2.010 

10.00 

upper  jump 

2.200 

9.00 

2.510 

7.30 

2.810 

6.00 

2.950 

/ 

3.140 

5.00 

3.460 

4.50 

3.770 

4.00 

4.090 

3.75 

. 

4.400 

3.50 

4.710 

3.25 

5.020 

/ 

4.75* 

v2nd  sub-harmonic 

5.150 

/ 

4.50 

upper  jump 

5.270 

/ 

3.50 

5.340 

2.50 

/ 
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(b)    Frequency  Decreasing: 


quency  of 

Amplitude  of 

Amplitude 

of   Remarks 

Input 

1st  Harmonic 

2nd  Sub-harmon 

tc 

(  c  ) 

( 

(    ) 

5 .  340 

2.50 

/ 

5.270 

/ 

3.50 

5.150 

/ 

4.50 

5.020 

/ 

4.75 

4.710 

/ 

5.50 

4.400 

/ 

6.50 

4.090 

/ 

7,80 

3.770 

/ 

9.50 

3.460 

/ 

11.00 

3.140 

/ 

13.50 

2.950 

/ 

14.50* 

*2nd  sub-hermonic 

2.810 

6.00 

lower  jump 

2.510 

7.25 

2.200 

8.75 

2.010 

10.00 

1.890 

11.00 

1.760 

12.25 

1.640 

13.00 

1.510 

15.00 

1 .  380 

17.00 

1.260 

21.20 

1.130 

24.00 

1.000 

25.50 

0.880 

0.45* 

*lst  Harmonic 

0.754 

0.33 

lower  jump 

0.628 

0,22 
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4-7  Analog  Computer  Analysis  for  a  3rd  Order  Sub -harmonic  Oscillation 
System: 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  A- 10: 


^itfy^U 


•  -r 


Fig.  4-10  Block  Diagram  of  a  3rd  Sub-harmonic 
Oscillation  System 

in  which  assume  the  nonlinearity  of  system  is  an  "ON-OFF"  element,  a 

sketch  of  the  characteristic  for  the  computer  simulation  is  shown  in 

Fig.  4-11,  it  is  approximated  from  three  straight  lines;  one  of  which 

with  a  slope  of  15.0;  one  of  which  with  a  slope  of  0.1;  and  the  other 

is  zero: 

The  circuit  setting  up  in  the  computer  is  shown  in  Fig.  4-12,  all 

components  to  be  used  are  schematically  indicated  with  the  diagram.   The 

characteristic  curve  from  the  computer  is  shown  in  Fig.  4-13. 

~-~'J*-T)    r — '  I — 


€L£ 


-H@- 


—4//..  ■ 


. 


/I 


! 
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Fig.  4-12  Function  Generator  Diagram  for  an 
"ON-OFF"  Nonlinearity 
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An  operational  diagram  cor  setting  up  the  analog  computer  is  shown 
in  Fig.  4-14,  in  which  the  value  of  all  components  to  be  used  are 
schematically  indicated  in  the  diagram.   The  scaling  factor  for  the 
computer  setup  is  the  same  as  diagram  of  Fig.  3-20  to  be  used. 

The  results  from  the  computer  are  condensed  in  Table  4-2,  and  with 
a  response  curve  shown  in  Fig.  4-15.   It  is  shown  that,  there  is  a  pair 
of  jumps  for  each  harmonic  oscillation  response  curve,  these  jump  pheno- 
menas  within  frequency  range  is  accompanied  shown  with  the  diagram  of 
Fig.  4-15. 
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Table  4-2:   Data  for  a  3rd  Order  Sub-harmonic  Oscillation  Response 

Curve  from  an  "ON-OFF"  Nonllnearity  System  and;  <.  =  0.4; 
K  ■   10.0;  9,-(t)  «  8.0  Cos(;ot  -»-  6) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Amp 

litdue  of 

Amplitude  of 

Remarks 

Input 

1st  Harmonic 

2nd 

Sub-har. 

3rd  Sub-har. 

(W  ) 

(  ft) 

(%) 

<r    ) 

0.628 

0.60 

0.880 

1.25 

1.130 

2.25 

1.380 

3.50 

1.640 

4.75 

1.760 

5.75 

1.890 

7.50 

2.010* 

50.00 

*lst  Harmonic 

2.140 

46.00 

Upper  Jump 

^.270 

42.50 

2.390 

39.00 

2.510 

35.00 

2.630 

31.50 

2.760 

28.00 

2.890 

26.00 

3.010 

22.50 

3.140 

20.00 

3.260 

/* 

*Transient  and 

4.270 

/* 

unstable 

4.390° 

"Transient  from 

4.520 

1st  Harmon* r 

4.650 

to  2nd  Sub-har. 

4.770 

4.990° 

5.020* 

5.140 

20.00 

*2nd  Sub-hai. 

5.270 

18.00 

Upper  Jump 

5.400 

16.70 

5.520 

16.00 

5.650 

15.00 

5.770 

13.80 

5.900 

13.50 

6.030 

12.70 

6.150 

12.30 

6.280 

11.90 

6.400 

11.30 

6.530 

11.00 

6.650 

10.60 

6.780 

10.20 

6.910* 

10.00 

*3rd  Sub-har. 

7.030 

35.00 

Upper  Jump 

7.160 

32.50 

7.280 

28.50 

7.410 

27.20 

7.530 

26.00 
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(b)      Frequency  Decreasing: 


Frequency  of 

Amplitude  of 

Amplitude  of 

Amplitude  of 

Pemarks 

Input 

1st  Harmonic 

2nd  Sub-har. 

3rd  Sub-har. 

(  ui>   ) 

(i  ) 

(  '  ) 

<   ) 

7.530 

26.00 

7.410 

27.20 

7.280 

28.50 

7.160 

30.50 

7.030 

32.50 

6.910 

35.00 

6.780 

36.50 

6.650 

37.50 

6.530 

40.00 

6.400 

42.50 

6.280 

43.75 

6.150 

47.00 

6.030 

50.00 

5.900 

52.50 

5.770 

56.00 

5.650* 

13.75 

*3rd  Sub-har. 

5.520 

15.00 

Lower  Jump 

5.400 

16.00 

5.270 

16.70 

5.140 

17.50 

5.020 

19.50 

4.900 

21.50 

4.770 

23.00 

.  4.650 

25.00 

4.520 

32.00 

4.390 

36.00 

4.27C  °* 

°*2nd  Sub-har. 

4.140 

Lower  Jump 

3.260° 

and  Transient 

3.140 

20.00 

Unstable  Region 

3.010 

22. 5< 

2.890 

26.00 

2.760 

27.50 

2.630 

31.00 

j 

2.510 

35.00 

2.390 

39.00 

2.270 

42.00 

2.140 

45.50 

2.010 

51.00 

1.890 

55.00 

1.760 

52.50 

1.640* 

5.00 

*lst  Harmonic 

1.5)0 

4.25 

Lower  Jump 

1.380 

3.50 

1.130 

2.30 

0.880 

1.25 

0.628 

0.60 
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CHAPTER  V 
FORCED  OSCILLATION  FOR  THE  HIGHER  ORDER 
FEEDBACK  CONTROL  SYSTEM 
5-1  General  Description: 

Up  to  the  Chapter  4,  only  forced  oscillations  in  the  case  of  the  2nd 
order  system  are  discussed.   In  the  extension  to  higher  order  systems , 
consider  a  block  diagram  shown  in  Fig.  5-1; 


Cy/V 


fie> 


Q(>) 


Fig.  5-1  Block  Diagram  of  a  Feedback 
Control  System 

in  which  G(s)  is  an  equivalent  transfer  function  of  the  linear  element 

of  the  system,  it  can  be  written  in  the  form: 


G(s)  «  K 


F(s) 


and  with  a  result  for  the  system: 

P(s)  E  +  Q(s)   Kf(£)  =  P(s) 

Where  T(s)  may  be  any  order  of  s;  as: 

~t  v  n       ,      n- 1       .      n-2  , 

P(s)  =  s     +  b.s         +  b-8  +.....+  b 

12  i 


_  ,   ,  m  m- 1  m-2 

Q(s)   =  s     +  c-S         +  c~s  + +  c 

1  l  m 


(5-1) 


(5-2) 


(5-3) 


Similarily  Q(s)   may  be 

I         m-9 

(5-4) 

In  general  the  frequency  responsive  element   2(s)  following  the  non- 
linearity  is  of  such  a  nature  as  to  attenuate  the  high  frequency,  that 
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means  the  order  of  P(s)  should  be  higher  than  the  order  of  Q(s).   The 
order  of  system  is  n.   Inserting  equations  (5-3)  and  (5-4)  into  (5-2)  and: 

(j&    f  h  i         .  )  (5-5) 

Equation  (5-5)  is  a  general  differential  equation  form  of  the 

forced  oscillation  of  any  order  feedback  control  system. 

5-2  Harmonic  Linearization  and  Iteration  Method  for  Investigating  the 
Equivalent  Gain  of  a  Nonlinearity : 

After  a  restoring  force  function  of  the  nonlinear  element  has  been 

defined,  the  equivalent  complex  gain,  N(E.)  of  the  nonlinear  element  is 

determined  by  the  ratio  of  the  complex  amplitude  of  the  first  harmonic 

of  the  output  to  the  amplitude  of  the  harmonically  varying  input  quantity. 

Recall  the  block  diagram  of  Fig.  5-1,  if  the  system  is  a  2nd  order 

3 
system  and  with  a  restoring  forcing  function;  f(E)  ■  a.E  -*-  a  E  ,  and: 

L10L1  -  «c  (5-6) 

s  (s  +«■•  )  ' 

or: 

E  +odE   +  K  N(E)  E  =  0,   +  U.  ®r  (5-7) 

Assume  the   input    function  is 

6>r  -  F  Cos(a)  t  +  9)  (5-8) 

and  one  solution  of  E: 

E  a  Ex  Cos  vd   t  (5=9) 

and  with  a  result: 

(  M)2  E  +  K  N(E1)E.)Cos^'t  -c/  E  Sin  .  t  ■  -A  Cos(^t  +  6) 

-  BSin(a  t  +  6)  (5-10) 

Comparing  the  result  with  the  equation  of  (3-42)  the  equivalent  gain: 

N(EX)  -  at  +|  a3E12  (5-11) 
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For  a  general  form  expression;  the  equivalent  gain  of  a  non- 
linearity  in  terms  of  the  amplitude  of  input  is  developed  by  Ya.  Z. 
Tsypkin.  The  result  for  a  symmetrical  nonlinearity  with  input  of  equation 
(5-9)  is: 

N(E1)  =  "if""  £  f(V   +  f(El/2)'  (5"12) 

If  more  accuracy  is  required,  the  equivalent  gain  N(E.)  may  be: 

N(E1)  s  3E1 1   £(E1)  +  f(Ei/2>  +  V3  f<   3  V2)j         <5-13) 

For  the  case  of  a  nonlinearity  with  a  restoring  force  function, 

3 
f(E)  ■  a  E  +  a  E  ,  the  results  from  equations  (5-12)  and  (5-13)  are  ex- 
actly the  same  as  equation  (5-11).  Hence,  for  any  nonlinearity  the  case 
of  first  harmonic  input  the  equivalent  gain  can  be  expressed  either  in  the 
form  of  equation  (5-12)  or  (5-13). 

Equation  (5-12)  or  (5-13)  not  only  substantially  simplifies  the 
calculation  for  a  second  order  forced  oscillation,  it  is  also  probably 
a  most  important  for  investigating  self  oscillation  or  forced  oscillations 
in  a  higher  order  system  with  single  or  multiple  nonlinearities  in  a 
feedback  control  system. 
5-3  Forced  Oscillations  for  a  Higher  Order  Feedback  Control  System: 


Consider  a  system  with  a  linear  transfer  function: 

K(s  +  cx) 
s(s  +  Pj)(s  +  P2)(s  +  p3) 


G(s)  ^         K(S  *  V  (5-14) 


or  in  the  form: 

6(.)  -   /<*  %C) 2 <5-15> 

s  4-   b  s  +  b2s  +  b  s 

From  the  equation  (5-5),  the  differential  equation  for  the  system 

of  equation  (5-15): 

(  — r  +  bj—y-  +  b2    2  +  b3  dr  dt 

dt       dt        dt 

4         3        2 

d?       V      2  F       3  ^  (5"16) 

01 


in  which; 


E  -  E.Cos    I 


■  FCos(   t  +  0) 

and; 

1      3 
f(E)  ■  &l   EjCos   t  +  7  a  E   (3Cos..t  +  Cos3 a  t) 

Inserting  these  equations  into  equation  (5-16)  and  vith  a  result: 
tO+-b,£,U)  t  Ka,  C£,  r  ^  K«-,  c  a  vc 

t    j  KQj  C  £* CQs3i*)&  -  ]p /ra3  £(J tv  A*«3 
~    Fw  [ vo  (a) *  -b*)&>S&  t  ( b,fij *"-  hj )./it*S>] 

(5-17) 
Neglect  the  higher  order  harmonic  terms  and  equate  the  coefficients 
of  Cos  !y  t  and  Sin  <  t: 

BCoO4-  -  b  f-jj?  Ms  c  *,  *J 

i  (U}*-b*) <*>&»&  •+  (b>i  14&0J  (5-18) 

—   ,.     tO*-b:.)iOA*ii&  -  CM        -C)  (5-19) 

Squaring  equations  (5-18)  and  (5-19)  and  adding: 

(5-20) 
Equation  (5-20)  is  the  equation  of  response  curve  of  the  system  of 
equation  (5-15).   If  the  constants  of  system  are  know,  a  response  curve 
in  the  £,  vfl)  plane  can  be  plotted. 

For  the  investigation  of  conditions  for  existence  of  forced  oscil- 
lation rearrange  the  equation  (5-20)  and  put: 
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3    „2 

al  +  4  *3E 


Hence 


<  / 


(5-21) 


or; 


t  [  U  •         K  (5-22) 

the  condition  for  existence  of  forced  oscillation  is  found  from 
equation  (5-22)  by  means  that  the  solution  of  to   should  be  real.   It  is 
to  be  noted  that  the  value  of  X  is  the  equivalent  gain  of  nonlinearity, 
its  value  should  be  larger  than  zero,  as  it  was  proved  in  Chapter  3. 

For  the  phase  investigating  recall  the  equation  (5-19)  and  re- 
write: 

•p 
==  .)  ~k  (<L,  (5-23) 


or: 


f I  ■  '  ~ t  (b 


hence 


in  which  the  value  of; 

=  Tim.  )]  (5- 

Equation  (5-25)  is  shown  that  the  value  of  6  is  a  function  of 
and  the  amplitude  of  first  harmonic  of  error  signal,  if  the  system 
constants  is  fixed. 
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If  we  use  the  equivalent  gain  for  the  system  Investigation,  and 
with  a  result: 

ft  —  Qr  (5-27) 

or : 


• 


I  r  Gcj'Mj/r  (a  (5-28) 

Where  E]  is  the  amplitude  of  first  harmonic  of  error  signal  of  forced 

oscillation. 

Equation  (5-28)  is  similar  to  a  linear  system,  the  techniques  for 

linear  system  are  applicable. 

5-4.   Analog  Computer  Analysis  for  a  Higher  Order  System: 
Consider  a  feedback  control  system  with  a  transfer  function  and  a 

nonlinearity  are  shown  in  Fig.  5-2: 


<g>-' 


7T" 

-ny 

± 


<Pc 


Fig.  5-2  Block  Diagram  for  a  Typical 
Feedback  Control  System 

in  which  the  restoring  force  function  is: 

f(E)  -  2E  +  0.020E3  (5-29) 

The  characteristic  curve  of  restoring  force  function  of  nonlinear- 
ity from  the  calculation  is  shown  in  Fig.  5-3.   Where: 

E  =  1.00;   2.00;   3.00;   4.00;   5.00;   6.00;   7.00;   7.50; 
f(E)  =  1.98;   3.84;   5.46;   6.72;   7.50;   7.68;   7.15;   6.60: 
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For  the  simulation  of  the  restoring  force  function  of  nonlinearity 
from  the  Computer,  it  is  approximated  by  three  straight  lines,  one  of 
which  with  a  slope  of  2;  one  of  which  with  a  slope  of  1;  and  the  other 
with  a  slope  of  zero.  The  circuit  setup  in  the  computer  is  shown  in 
Fig.  5-4,  all  components  used  are  schematically  in  the  diagram. 

The  characteristic  curve  of  restoring  force  function  from  the  com- 
puter is  shown  in  Fig,  5-5;  (a)  shows  the  waveforms  of  the  input  and  the 
output;  (b)  shows  the  characteristic  curve  of  the  nonlinearity. 


lit*)  / — \   Ekv 

-  :  ■ 


(fy — 


T  £** 


Fig,  5-4  Analog  Computer  Setup  for 

f(E)  =2E  ■  0.020E3  Simulation 

The  operational  block  diagram  and  the  operational  circuit  diagram 
are  shown  in  Fig.  5-6  and  Fig.  5-7.  The  all  components  used  in  the 
operational  circuit  diagram  are  sechematically  in  the  diagram. 

The  results  from  the  computer  are  condensed  in  Table  5-1,  and  with 
response  curves  shown  in  the  £v0*  plane  and  fiL'^>  plane  are  shown  in  Fig.  5-8 
and  Fig.  5-9  respectively. 
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Fig.  5-6  Operational  Block  Diagram  setup 
in  the  Analog  Computer 
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Table  5-1.   Data  for  Response  Curve  from  Computer  for  a  Higher  Order 

System;  in  which  f(E)  -  2E  -  0.020E3,  G,  (t)  -  3.0  Cosfc't  +G) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(M   ) 

(  £?/) 

(  a  ) 

( 

( 

0.314 

0.40 

270.00 

3.10 

0.00 

0.377 

0.55 

270.00 

3.30 

0.00 

0.440 

0.75 

262.00 

3.40 

0.00 

0.503 

0.90 

255.00 

3.50 

5.80 

0.565 

1.10 

246.00 

3.60 

9.00 

0.628 

1.25 

241.00 

3.80 

12.50 

0.755 

].70 

235.00 

4.30 

20.00 

0.880 

2.50 

220.00 

5.00 

24.50 

1.000 

7.40 

105.00 

6.30 

107.50 

1.130 

6.90 

64.00 

5.10 

120.00  ' 

1.260 

6.50 

37.50 

4.25 

137.50 

1.380 

6.10 

25.00 

3.50 

155.00 

1.510 

5.60 

15.00 

3.00 

165.00 

1.640 

5.25 

10.00 

2.50 

172.00 

1.760 

4.90 

7.50 

2.00 

180.00 

1.890 

4.50 

5.00 

1.70 

186.00 

2.200 

4.00 

0.00 

1.20 

195.00 

2.500 

3.50 

0.00 

0.85 

210.00 

2.810 

3.30 

0.00 

0.60 

220.00 

3.140 

3.20 

0.00 

0.40 

232.00 

3.460 

3.10 

0.00 

0.30 

245.00 

3.770 

3.00 

0.00 

0.25 

260.00 

4.090 

3.00 

0.00 

0.20 

270.00 

(b)  Frequency  Decreasing 


Frequency  jf 

Amp 

litude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(a.  ) 

(  *,) 

(  ':  ) 

(  -  ) 

(  6  ) 

2.510 

3.50 

0.00 

0.78 

210.00 

2.200 

4.00 

0.00 

1.20 

195.00 

1.890 

4.50 

5.00 

1.65 

185.00 

1.760 

4.90 

7.50 

2.00 

180.00 

1.640 

5.25 

10.00 

2.50 

171.00 

1.510 

5.70 

15.00 

3.00 

162,00 

1.380 

6.10 

25.00 

3.50 

155.00 

1.260 

6.60 

37.50 

4.25 

137.00 

1.130 

6.60 

64.00 

5.00 

121.00 

1.000 

7.40 

105.00 

6.25 
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CHAPTER  VI 
CONCLUSIONS 

The  aim  of  this  dissertation  was  to  present  a  method  for  investigat- 
ing the  conditions  for  existence  of  the  forced  oscillations  of  a  nonlinear 
feedback  control  system.   This  method  is  called  "harmonic  linearization 
method",  it  is  based  on  the  harmonic  balance  and  iteration  method.   From 
this  method,  the  conditions  for  existence  of  the  fundamental  frequency 
and  the  sub-harmonic  of  order  2  and  3  are  investigated,  and  also  a  equi- 
valent gain  of  the  nonlinearity  in  terms  of  input  amplitude  and  frequency 
has  been  developed. 

The  term  of  the  "response  curve"  as  used  here  is  defined  as  the  combina- 
tion of  two  components,  one  is  the  attenuation  frequency;  it  is  either 
shown  in  the  "Error  vs  Frequency"  plane  or  "Output  vs  Frequency"  plane, 
and  the  other  is  the  phase  frequency  response,  it  is  with  a  difference 
value  from  the  different  representation  of  attenuation  frequency  response. 

The  results  of  this  analysis  show  that  a  nonlinear  feedback  control 
system  may  possess  more  than  one  type  of  forced  oscillations,  which  type 
will  exist,  it  is  depending  on  the  characteristic  of  the  system  and  the 
characteristic  of  input  forcing  function. 

For  any  nonlinear  feedback  control  system,  the  response  curve  of 
either  fundamental  or  sub-harmonic  forced  oscillations,  it  is  usually  a 
sudden  jump  in  amplitude  and  phase  when  a  jump  transient  take  place.  The 
range  of  jump  frequency  depends  on  the  value  of  damping  of  the  system  and 
the  amplitude  of  the  input;  usually  it  is  increased  as  the  damping  decreas- 
ing, and  decreased  as  the  input  amplitude  decreasing.  A  vertical  tangents 
locus  investigating  for  calculating  the  jump  frequencies  in  terms  of  system 
constants  and  the  amplitude  of  input  was  investigated. 
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Due  to  the  resonance  phenomena  and  the  linear  approximations  of  re- 
storing force  function,  there  is  some  difference  between  the  theoretical 
and  the  experimental  response  curves.   It  If  only  exlstant  in  the  jump 
frequency  range,  other  than  this  range  of  frequencies,  it  is  very  close 
to  each  other. 

Sub-harmonic  oscillations  are  undesired  for  a  nonlinear  feedback 
control  system,  it  is  actually  an  unstable  phenomenon.   Under  suitable 
conditions,  a  multiple  order  sub-harmonic  forced  oscillation  may  exist 
in  a  system  with  a  constant  input  forcing  function.   There  is  also  a 
jump  phenomena  for  each  order  of  sub-harmonic  oscillations.   The  condi- 
tions for  the  existence  of  the  2nd  and  3rd  order  sub-harmonic  forced 
oscillations  are  investigated,  and  also  theoretical  response  curve  equa- 
tions for  each  case  are  developed. 

It  is  possible  an  important  advantage  that  the  design  and  analysis 
of  a  nonlinear  feedback  control  system  by  using  the  harmonic  lineariza- 
tion method,  by  means  for  investigating  an  equivalent  gain  of  the  non- 
linearity.   In  this  case,  all  principles  and  techniques  applied  for  linear 
system  are  applicable. 
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